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ABSTRACT 

We have obtained low resolution spectra of nineteen red and blue low surface brightness galaxies, 
using the Marcario Low Resolution Spectrograph on the 9.2m Hobby-Eberly Telescope. These galaxies 
form a very heterogeneous class, whose spectra qualitatively resemble those of high surface brightness 
galaxies covering the full range of spectra seen in galaxies of Hubble types from E to Irr. We use a 
combination of emission line (EW(Ha), [N II]A6584/Hq;) and absorption line (Mgb, H/3, <Fe>) based 
diagnostics to investigate the star-formation and chemical enrichment histories of these galaxies. These 
are diverse, with some galaxies having low metallicity and very young mean stellar ages, and other 
galaxies showing old, super-solar metallicity stellar populations. In contrast with some previous studies 
which found a strong trend of decreasing metallicity with decreasing central surface brightness, we find 
a population of galaxies with low surface brightness and near-solar metallicity. Correlations between 
several of the gas phase and stellar population age and metallicity indicators are used to place contraints 
on plausible evolutionary scenarios for LSB galaxies. The rcdshift range spanned by these galaxies is 
broad, with radial velocities from 3400 km/s to more than 65000 km/s. A subset of the sample galaxies 
have published H I redshifts and gas masses based on observations with the Arecibo 305m single-dish 
radio telescope, which place these galaxies far off of the mean TuUy-Fisher relation. Our new optical 
redshifts do not agree with the published H I redshifts for these galaxies. Most of the discrepancies can 
be explained by beam confusion in the Arecibo observations, causing erroneous H I detections for some 
of the galaxies. 

Subject headings: galaxies: abundances — galaxies: distances and redshifts — galaxies: stellar content 
— galaxies: fundamental parameters — galaxies: formation and evolution 



1. INTRODUCTION 

The class of low surface brightness (LSB) galaxies is 
a diverse one. LSB galaxies exhibit a range of absolute 
magnitudes, colors, and morphologies comparable to those 
of high surface brightness (HSB) galaxies. The defining 
characteristic of LSB galaxies is, of course, their faint cen- 
tral surface brightness, variously defined as being at least 
one magnitude fainter than the cannonical Freeman (1970) 
value of 21.65 B mag/D", or fainter than 23 B mag/D"; 
the precise limit varies from investigator to investigator. 
However, whether the low surface brightness results from 
a common evolutionary history for these galaxies remains 
unknown. 

Studies of LSB galaxies fall into two broad classes. The 
first class comprises the surveys to find LSB galaxies, and 
attempts to quantify what fraction of the total galaxy 
number density, luminosity density, and mass density are 
made up by them. The second class comprises studies of 
the internal structure of LSB galaxies, including investi- 
gations into their gas content, star formation history, and 
evolutionary processes. 



The large area surveys which were initially mined for 
LSB galaxies were photographic surveys, including the 
Upsala Galaxy Catalog (UGC; Nilson 1973), the POSS- 
2 survey (Schombert & Bothun 1988; Schombert et al. 
1992), and the APM survey (Impey et al. 1996). In the 
past decade, large area CCD surveys have also found a 
significant population of LSB galaxies. Notable among 
these are the drift scan survey of Dalcanton et al. (1997), 
who use well quantified automated detection methods, and 
the multicolor survey of O'Neil, Bothun, & Cornell (1997, 
hereafter OBC97) and O'Neil et al. (1997b, hereafter OB- 
SCI97), which turned up a previously unseen population 
of red LSB galaxies. 

These surveys proved the existence of LSB galaxies as 
more than the occasional oddity. Dalcanton et al. (1997) 
make a conservative estimate of the number density of LSB 
galaxies of A = O.Olto^oos galaxies /iggMpc^^, equal to the 
number density of HSB galaxies. However, the Dalcanton 
et al. survey suffers from small number statistics, with 
only 7 LSBs detected. McGaugh, Bothun, & Schombert 
(1995) have a larger sample, and come to a similar conclu- 
sion about the LSB galaxy number density. More recently. 



1 



2 



de Jong & Lacey (2000) have used a photographic survey 
of ~ 1000 Sb-Sdm galaxies to calculate the bivariate lu- 
minosity and number density of galaxies as a function of 
both surface brightness and scale length. They find that 
while the luminosity density is dominated by HSB galax- 
ies (peaked at /xb(0) ~ 21.0 mag/D"), the number density 
of galaxies has an essentially flat distribution at all sur- 
face brightnesses fainter than fisiO) — 21.0 mag/D", and 
therefore LSB galaxies make up a significant fraction of the 
total number of galaxies (See also the surveys of O'Neil & 
Bothun (2000) and Cross & Driver (2002)).' 

The best studied LSB disk galaxies are those in the the- 
sis sample of de Blok (1997), originally selected from the 
photographic survey catalogs of Schombcrt et al. (1992) 
and the UGC catalog. Observations of these ^ 20 galax- 
ies include multicolor optical surface photometry (de Blok, 
van der Hulst, & Bothun 1995), VLA mapping of the H I 
distribution and kinematics (de Blok, McGaugh, & van der 
Hulst 1996; de Blok & McGaugh 1997), and spectroscopy 
of the H II regions (de Blok & van der Hulst 1998). The 
general picture obtained from these observations is that 
LSB galaxies are metal poor and gas rich, with a range 
of sizes (scale lengths) comparable to HSB galaxies, but 
with much lower surface densities. However, note that de 
Blok's selection criteria included an existing H I detection. 
The kinematics, and therefore the gravitational potentials, 
are dominated by dark matter at all radii (de Blok & Mc- 
Gaugh 1997; McGaugh, Rubin, & de Blok 2001), and the 
low global star formation rates (SFR) are the combined 
result of both the low surface density, and the low metal- 
licity (Gerritsen & de Blok 1999). 

Sprayberry, Impey, Bothun, & Irwin (1995) obtained 
photometry, low-resolution optical spectroscopy, and HI 
spectroscopy for a sample of 8 giant LSB galaxies. Their 
optical spectroscopic data pertain primarily to the bulges 
of these galaxies, and they find no difference between the 
bulges of large LSB and HSB galaxies. 

The wide field CCD survey of OBC97 found a large num- 
ber of LSB galaxies, covering a wider range of colors than 
those found by previous photographic surveys. In partic- 
ular, red LSB galaxies with (B— V)> 0.9, colors consistent 
with purely old stellar populations, were catalogued for the 
first time. The follow-up single dish H I observations us- 
ing the Arecibo 305m single-dish radio telescope (O'Neil, 
Bothun, & Schombert 2000a) suggested that many of these 
red LSBs (and some of the blue LSBs) were quite pecidiar 
indeed, having gas mass-to-light ratios (M/L) in excess 
of 8, with a few higher than 30. Further, the rotational 
line widths were far larger than would be expected from 
their luminosities given the standard TuUy-Fisher (1977) 
relation for spiral galaxies. 

Larger scatter in the TuUy-Fisher relation among LSB 
galaxies was also measured by Burkholder, Impey, & 
Sprayberry (2001). Their survey of over 250 HSB and LSB 
galaxies selected from the APM survey (Impey ct al. 1996) 
included photometry, HI spectroscopy, and low resolution 
optical spectroscopy (Impey, Burkholder, & Sprayberry 
2001; Burkholder, Impey, & Sprayberry 2001). They used 
emission line ratios to measure the metallicity for 49 HSB 
and 17 LSB galaxies (where they define /Ub = 22mag/n"as 
the dividing line) and found there to be a roughly 50% 
overlap in the metallicity range of the two groups. 



Bell et al. (2000) combined optical and near-infrared sur- 
face photometry to study the stellar populations in subsets 
of both the de Blok sample and the red galaxies from the 
OBSCI97 survey. Broadband optical studies of stellar pop- 
ulations suffer from a degeneracy in color space between 
age and metallicity; red colors could indicate either older 
stars or more metal-rich stars, while bluer colors indicate 
either younger stars or metal-poor stars. The inclusion of 
near-infrared colors partially breaks this degeneracy, es- 
sentially giving a handle on the fractional light contribu- 
tion from old stars which dominate the light in the near- 
infrared, and younger stars which will dominate at bluer 
wavelengths. Bell et al. found that the bluer galaxies from 
the de Blok sample were younger than similar HSB late- 
type galaxies. Further, the five red LSB galaxies in their 
study were a mixture of two types: genuinely red galaxies 
with old stellar populations, but central surface bright- 
nesses on the bright end of the LSB range 22.5mag/n" 
in B), and a group of objects with lower surface brightness 
whose true colors, after correcting for galactic extinction, 
were in fact as blue as the de Blok sample galaxies and 
had similar stellar mean ages. Interestingly, the three red 
galaxies with old stellar populations were not detected in 
the O'Neil et al. (2000a) H I survey, while the two galaxies 
with blue extinction corrected colors were detected. How- 
ever, none of the red galaxies with extreme Mhi/L ratios 
were studied by Bell et al. 

We have obtained spectra for a sample of 19 red and 
blue LSB galaxies chosen from the catalog of OBC97, with 
the intent of studying the star formation histories of these 
galaxies. Using a combination of emission and absorption 
line measures, we compare the chemical enrichment of the 
stars and gas, and we use the stellar population models 
of Vazdekis (1999, 2000^, hereafter V1999 and V2000, re- 
spectively) to interpret these data quantitatively. In par- 
ticular, we attempt to determine whether a single LSB 
galaxy formation model can be reconciled with the vari- 
ous observational attributes of these galaxies. 

This paper is organized as follows. In Section 2 we 
present the sample, the spectroscopic and photometric ob- 
servations, and the data reduction. Sections 3 through 
8 presents the LSB spectra, analysis of the emission line 
and absorption index strengths, and comparison to HSB 
galaxies. We discuss these results in the context of LSB 
formation models in Section 9. Finally, we summarize our 
conclusions in Section 10. 

2. OBSERVATIONS & DATA REDUCTION 

2.1. Sample selection for the spectroscopic observations 

This project began during the early science operations of 
the Hobby-Eberly Telescope (HET; Ramsey et al. 1998) in 
1999 November and observations continued through 2001 
January. For the first set of spectroscopic observations, 
during winter 1999-2000, we chose a group of LSB galaxies 
from the sample of OBC97 which covered a range of color 
(0.3 < B— V < 1.35) and central surface brightness (22.0 
< /u_b(0) < 24.0). Twelve galaxies were observed during 
the first observing period. A new sample of six galaxies 
was chosen for observation during winter 2000-2001. These 
objects were selected from the H I follow-up survey of 
O'Neil et al. (2000a), and stand out due to their extremely 

1 http; / /www.iac.es/galeria/vazdekis/MODELS_2000/out_li_BI 
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Table 1 



LiTERATrRE DATA FOR CALAXIES IK SPECTROSCOPK' SAMPLE" 



galaxy 


RA 


Dec 


cz (HI) 




B(5 ) 


L.(5 ) 


Btot 


T 1' 
ictot 


b 

Mb(0)cui- 


■ c 

I 


E(B-V) 






(J2000) 


(J2000) 


(km s ) 


(") 


(mag) 


(mag) 


(mag) 


(mag) 


(mag/D") 


(dog) 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 




OBC971P1-3 


■ ^ i- ■ i-O 


Uo. -LI. 


3746 


7 1 


18 28 


16 70 


17 15 


15 81 


22 2 


45 


071 




OBC971p6-1 




uO.O 1 . 


10882 


18.9 


19.55 


16.79 


17.23 


15.04 


24.2 


66 


0.061 




OBC97 P9-4 


OQ.1 Q.AC\ 




4205 


10 9 


20 84 


18 23 


18 86 


16 29 


25 


67 


144 




OBC97 Cl-2 


<JO . l-iJ .'JO 


, r_l -L . 


8531 


5.2 


18.99 


16 91 


18 10 


15 84 


22 5 


21 


049 


[OBC97]Cl-4 




91 -nn- 1 9 




12 6 


18 09 


15 80 


16 55 


14 48 


22 3 


63 


0-50 


[OBC97]C3-2 




20:59:40 




6.0^^ 


18.17 


1d.uo±u.0d 


16.34 


14. 14 


22.8 


47 


0.038 


pBC97 C4-1 


08:24:33 


21:27:04 


7905 


11. 0'^ 


18.97 




17.83 




23.6 


39 


0.044 




OBC97 


C4-5 


08:25:48 


21:52:45 




5.6 


18.96 




18.23 




23.0 


57 


0.046 




OBC97 


C5-3 


08:27:44 


22:28:52 


12942 


8.6 


20.12 


18.31±0.1 


18.88 


16.37 


24.0 


20 


0.041 




OBC97 


C5-5 


08:27:12 


22:53:39 


5524 


23.0 


19.90±0.2 


18.60±0.2 


18.08±0.2 


16.19±0.2 


24.3 


45 


0.038 




OBC97 


C6-1 


08:27:31 


21:30:03 


4322 


6.5 


18.57 


16.22±0.06 


17.53 


15.27 


22.5 


44 


0.041 




OBC97 


C8-3 


08:17:17 


21:39:12 


3570 


4.8 


19.22 




17.69 




24.3 


72 


0.039 




OBC97 


N3-1 


12:32:59 


07:29:26 




8.6 


19.05 




18.02 




23.5 


64 


0.023 




OBC97 


N9-1 


10:20:18 


27:45:10 




7.6 


18.41 


19.18±1 


17.29 


17.60±1 


22.3 


41 


0.036 




OBC97 


N9-2 


10:20:22 


28:07:56 


7746 


6.4 


19.80 


20.59±0.06 


18.89 


17.49 


24.0 


47 


0.037 




OBC97 


NlO-2 


11:58:42 


20:34:43 




10.8 


18.03 


16.22±0.06 


16.57 


14.85 


22.8 


52 


0.023 




OBC97 


NlO-4 


11:58:52 


20:58:28 


7478 


9.6 


19.94 




18.90 




24.2 


31 


0.027 




OBC97 


Ul-4 


11:38:25 


17:05:14 


3450 


32.5 


19.00 


17.26±0.06 


13.08 


12.21 


23.3 


55 


0.025 




OBC97 


Ul-8 


11:39:23 


16:50:31 




10.3 


19.30 


16.55±0.06 


18.15 


15.65 


24.1 


70 


0.025 



'^Columns (l)-(3) give the names and positions from O'Ncil et al. (2000a), and the rcdshifts of their HI detections are listed in column 
(4). The exponential scale lengths measured by OBC97 have been converted to half light radii using the formula r^ — 1.68a (Sparks 
1988), and these are listed in column (5). Columns (6)-(7) list the B-band and I-band magnitudes inside apertures of radius 5", and 
the corresponding total magnitudes are given in columns (8)-(9) (OBC97, OBSCI97). The inclination corrected B band central surface 
brightnesses are given in column (10), and the inclinations are in column (11), both from OBC97. Column (12) lists the galactic 
extinction derived from the Schlegel et al. (1998) dust maps. 

^Unless otherwise noted, the typical intrinsic uncertainty on the photometry is reported to be ±0.05 mag for the aperture magnitudes, 
and ±0.1 mag for the total magnitudes and /^B(0)cor. The uncertainty in the photometric zeropoints are 0.03 mag in the B-band, and 
0.04 mag in the I-band (OBC97, OBSCI97). 

*^The intrinsic uncertainty on the inclination angle is ±5° (OBC97). 

'^Not measured by OBC97. We use this value for spectral extraction to optimize the signal-to-noise ratio, (see Sect. 2.3.3). 



high Mhi/Lb, and that they he more than 5 a from the 
I-band TuUy-Fisher relation (Tully & Fisher 1977; Pierce 
& Tully 1988). Additionally, one galaxy ([OBC97]P06-l2) 
was chosen because it is the first LSB (albeit a giant LSB) 
galaxy for which CO has been detected (O'Neil, Hofner, & 
Schinnerer 2000b; see also Matthews & Gao 2001). These 
seven galaxies fall into the same range of color and central 
surface brightness as the original sample. We obtained 
spectra for a total of 19 galaxies. Table 1 lists the galax- 
ies in our spectroscopic sample, along with positions and 
relevant data from the literature. 

2.2. Photometry 

2.2.1. Observations & Basie Reductions 

We obtained new photometry for a subset of the LSB 
galaxies in the spectroscopic sample, using the Prime Fo- 
cus Camera (PFC; Claver 1995) on the McDonald Ob- 
servatory 0.8m telescope in November and December of 
1999. The PFC is a wide-field imager, with a 45' x 45' 
field of view. The wide-field images include eleven of the 
nineteen galaxies in the spectroscopic sample, as well as 
sixteen galaxies from the OBC97 catalog not included in 
our spectroscopic observations. The instrumental config- 
uration is given in Table 2. Images were taken through 
the B, Rc, and Ic filters, centered on the OBC97 survey 
fields CI, C3, C4, C5, N9, and Ul. Multiple images were 
taken in each band to aid in cosmic ray rejection. Ad- 

^ The official (lAU) names for the LSB galaxies in this work origi- 
nated in OBC97, and are listed in full in Table 1. For convenience, 
however, we will use abbreviated galaxy names without [OBC97] for 
the remainder of this paper. 



ditionally, we took images of 10 fields containing a total 
of 217 standard stars included in Landolt (1992) to facil- 
itate photometric calibration of the data. For each field 
observed. Table 3 lists the position of the field center, the 
total integration time and number of exposures through 
each filter, and the image quality of the coadded data, de- 
scribed in terms of the FWHM and ellipticity of the PSF 
as determined from gaussian fits to radial profiles of more 
than 100 stars per image. 

Data reduction was done within the IRAF^ environ- 
ment. The images were overscan and bias subtracted, 
then a shutter correction was applied. Next, the images 
were divided by a normalized dome fiatfield, as well as a 
normalized illumination correction. The illumination cor- 
rection is a highly clipped average of all science frames 
observed at low airmass, in the absence of clouds or moon- 
light, and corrects large scale gradients not reproduced by 
the dome flatfield image. After this step, the science im- 
ages were individually examined for residual real gradients 
caused by scattered moonlight, clouds, or sky illumination 
(see Chromey & Hasselbacher 1996). Any residual gra- 
dients larger than 1% across the full image were fit using 
imsurf it and subtracted off, with the mean sky level then 
added back in as a constant to maintain counting statis- 
tics. 

Prior to coaddition, a badpixel mask was generated for 
each image by taking the standard mask file for the CCD, 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 2 
Instrumentation 



Spectroscopy Photometry 



Telescope 


Hobby-Eberly Telescope (9.2m) 


McDonald Observatory 0.8m 


Instrument 


Low Resolution Spectrograph 


Prime Focus Camera 


Dates 


10/99-2/01 


11/99-12/99 


CCD 


Ford Aerospace 3072 X 1024 


Loral-Fairchild 2048 X 2048 




(binned 2x2) 




Gain 


1.82c" /ADU 


1.6 e-/ADU 


RON 


7.0 e" 


5.87e- 


Pixel size spatial 


0.465" 


1.354" 


Wavelength Range 


4300A - 7240 A 


B, Rc, Ic 


Slit Width X Length 


2"x 4' 




Grism 


600 1/mm 




Dipi")orsion 


2A/i)ix 




lirSoluLiull(cT) 


:i.:!A (l!Jl) km ' <! 52l)()A) 





Table 3 

PFC Photometry Observations 



Field* 


Filter 


Exposure Time*^ 


Nobs" 


FWHMpsf'^ 


epSF° 


NpSF* 


CI 


B 


300 


1 


2.8 


0.12 


176 


(08:19:52.6 +21:00:17) 


Rc 


180 


1 


2.7 


0.13 


111 




Ic 


120 


1 


2.4 


0.12 


101 


C3 


B 


1500 


3 


3.0 


0.08 


229 


(08:22:46.6 +21:03:34) 


Rc 


780 


3 


2.7 


0.06 


176 


Ic 


600 


3 


3.0 


0.04 


127 


C4 


B 


2100 


4 


2.8 


0.05 


508 


(08:26:01.3 +21:57:06) 


Rc 


1380 


5 


2.8 


0.04 


201 


Ic 


600 


3 


3.0 


0.06 


109 


C5 


B 


2100 


4 


3.1 


0.10 


180 


(08:26:54.1 +22:35:23) 


Rc 


1080 


4 


3.0 


0.04 


126 


Ic 


600 


3 


3.1 


0.07 


145 


N9 


B 


2100 


4 


3.2 


0.07 


264 


(10:20:29.3 +27:55:29) 


Rc 


780 


3 


2.8 


0.07 


141 




Ic 


600 


3 


3.0 


0.05 


118 


Ul 


B 


1200 


2 


3.8 


0.04 


147 


(11:39:29.1 +17:04:21) 


Rc 


1200 


4 


3.2 


0.03 


101 


Ic 


720 


3 


3.5 


0.06 


145 



^First row of each entry is the field name from OBC97. Second row lists the coordinates (J2000) 
of the field center. 

^Total exposure time for all observations, in seconds. 

"^Number of individual exposures. 

^Mean full-width at half-maximum of a gaussian fit to the point-spread-function (PSF), measured 
in arcsec. 

^Mean measured ellipticity of the PSF. 

^Number of stars used in determination of FWHM and e. 
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and manually adding all the pixels affected by satellite 
trails or by ghost images caused by reflections within the 
PFC optics. The coaddition process was done using the 
task imcoadd (J0rgensen, in prep.) in the Gemini IRAF 
package. This task determines spatial offsets between all 
the input images, shifts them accordingly, generates cos- 
mic ray masks, and averages all the good pixel values to 
produce a coadded output image. The photometric zero- 
point offsets between all the input images and the output 
image; were determined. As long as at least one input im- 
age was taken during photometric conditions, the coadded 
output can then be easily calibrated onto a standard sys- 
tem. Finally, we used the HST guide star catalog reference 
frame to establish the world coordinate system (WCS) for 
the coadded images. 

2.2.2. Object Extraction & Standard Photometry 

Object extraction was done using SExtractor (version 
2.1.6; Bertin & Arnouts 1996), and the output catalogs 
were matched using the WCS positions to produce a single 
table with photometry in all three bands (B, Rg, Ic) for 
each object. 

To calibrate the photometry onto a standard system, we 
used observations of the Landolt (1992) selected area fields 
SA92, SA95, SA98, SAlOl, SA104, and SA113, as well as 
the fields near the standards PG0918-h029, PG0231-h051, 
PG1633-I-091, and PG0233-I-051. Observations were made 
each photometric night, interspersing science and standard 
fields. In total, we have 380 stellar observations in each of 
the three bands, covering 217 different stars. The standard 
fields were observed over a range of airmasses spanning 
from 1 to 3, to better constrain the extinction coefficients. 
In contrast, the science fields were all observed in the air- 
mass range 1 to 1.25. Average extinction coefficients cover- 
ing all three photometric nights were used in the standard 
calibration. The rms scatter of the final standard calibra- 
tion of the photometric zeropoints was 0.06 mag (B), 0.05 
mag (Rc), and 0.05 mag(Ic). 

There were no observations of field Ul taken during 
photometric conditions. To calibrate this field, we used 
the ^ 1000 stars in our observation in common with the 
USNO-A2.0 catalog (Monet et al. 1998). We derived rela- 
tions between the SExtractor "best" magnitudes and the 
USNO magnitudes for both the Ul and N9 fields, and 
used the zeropoint offset between these two relations to 
calibrate the photometry for the Ul field to the standard 
system. The USNO magnitudes do not have to be on a 
standard system to make this technique work, merely be 
internally consistent. The extra uncertainty contributed 
by this calibration is 0.02 mag, smaller than our zeropoint 
uncertainty as derived using SExtractor. 

We applied galactic extinction corrections to our pho- 
tometry based on reddening values derived from the 
Schlegel et al. (1998) dust maps, and the coefficients de- 
rived by them for extinction in each band versus color 
excess (e.g. Ab/E(B-V) = 4.315, Av/E(B-V) = 3.315, 
Ar^/E(B-V) = 2.673, AijE(B-V) = 1.940). Table 4 
lists the extinction corrected magnitudes for the 27 galax- 
ies in the OBC97 sample that arc included in our obser- 
vations with the PFC. Also listed are the E(B-V) values 
used to determine the extinction correction in each band. 

2.2.3. Comparison to literature data 



We have compared our photometry with the measure- 
ments of OBC97 and OBSCI97 (which were also obtained 
with the PFC). Figure 1 shows the comparison. 

Figures la and Ic show the comparison of the aperture 
magnitudes in the B-band and I-band, respectively. In 
the B-band, we find an offset of —0.52 ± 0.04 mag with an 
rms scatter of 0.19 mag (excluding the two most deviant 
measurements, N9-2 and C4-2). Offsets are our magni- 
tudes minus OBC97 magnitudes. In the I-band, we find 
an offset of —0.31 ± 0.08 mag with an rms scatter of 0.34 
mag (excluding N9-1 and N9-2). The offsets and the large 
scatter may be due to the sensitivity of the small apertures 
used (5") to changes in seeing or image quality. In the in- 
terim between when OBC97 obtained their data and the 
time when we made our observations, the PFC optics were 
anti-reflection coated and re-aligned, which improved the 
delivered image quality. In our coadded images, the stel- 
lar PSF has a typical FHWM of 2.2 pixels (2.9"), with an 
ellipticity of 0.07. O'Neil has provided some of her images 
to the NED^ database, and the typical image quality we 
measure for these small images is ~ 3 pixels fwhm (4"). 
Poorer image quality will scatter more of the light outside 
the small 5" aperture, resulting in fainter aperture magni- 
tude measurements, consistent with the offsets we see. 

OBC97 derive total magnitudes from exponential and 
King proflle fits to the galaxy surface brightness profiles. 
We compare these to the SExtractor "best" magnitudes, 
which arc also a measiirc of the total magnitude (though 
systematically too faint by 6%, see Bertin & Arnouts 
(1996) for details) in Figures lb and Id. The agreement in 
the B-band is good, except for two outliers, N9-2 and Ul-4. 
Examination of the surface brightness profiles in OBC97 
for these two objects shows the profiles to be more complex 
than could be well fit by a single exponential. Excluding 
these two points, the offset between the two datasets is 
—0.11 ± 0.04 mag with an rms scatter of 0.21 mag. In the 
I-band, the offset is 0.23 ± 0.08 mag with an rms scatter 
of 0.33 mag. 

We have not measured central surface brightness, /x(0), 
from our photometry, so we cannot compare to the values 
in OBC97. However, since they derive total magnitudes 
based on /i(0) and the exponential scale length (a, derived 
over an area with radius typically larger than 12"), we 
expect offsets in /x(0) to be comparable to offsets in the 
total magnitudes, which are small in the B-band. 

Figures Ic If show the comparisons of the (B— Ic) aper- 
ture and total colors, respectively, with the results from 
OBSCI97. For both the aperture and total colors, we de- 
rive a zeropoint offset of —0.24 ± 0.07 mag with an rms 
scatter of 0.30 mag, when we exclude the two most de- 
viant points (N9-1 and N9-2 for the aperture colors, N9-1 
and C5-3 for the total colors). For N9-1 the uncertainty on 
(B— Ic) from OBSCI97 is 1 mag. The typical intrinsic un- 
certainty for our (B— Ic) color measurements is 0.05 mag, 
derived from the SExtractor uncertainties on the individ- 
ual band magnitudes, added in quadrature. In the follow- 
ing analysis we use our new photometry when available. 
For the remaining galaxies with spectroscopic observations 
we use photometry from OBC97. We calibrate the colors 
from OBSCI97 to be consistent with our data by applying 
an offset of —0.24 mag, and adopt an uncertainty in this 

^ http://ned.ipac.caltech.edu 
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Table 4 

New Photometry from PFC Obser\'atio\s" 



Name 



B(5") 



Ro(5" 



Io(5" 



Bt. 



lie. tot 



Ic.tot 



E(B-V) 



Cl-1 
Cl-2 
Cl-4 
Cl-6 
C3-1 
C3-2 
C3-6 
C4-1 
C4-2 
C4-3 
C4-5 
C4-6 
C4-7 
C4-8 
C5-1 
C5-2 
C5-3 
C5-4 
C5-5 
N9-1 
N9-2 
Ul-1 
Ul-2 
Ul-3' 
Ul-4' 
Ul-6' 
Ul-8' 



17.95±0.03 
18.66±0.06 
17.45±0.02 
17.53±0.02 
18.26±0.01 
17.60±0.01 
17.81±0.01 
18.58±0.03 
19.05±0.05 
18.46±0.03 
18.62±0.03 
17.65±0.01 
19.63±0.06 
1!).11±0.05 
17.()5±0.01 
18.91±0.02 
19.63±0.04 
18.14±0.01 
19.28±0.03 
18.01±0.01 
20.32±0.07 
18.17±0.02 
19.19±0.04 
17.96±0.01 
18.05±0.02 
19.32±0.04 
18.73±0.02 



16.31±0.01 
17.26±0.03 
16.01±0.01 

16.48±0.01 
16.63±0.01 
16.24±0.00 
16.60±0.01 
18.00±0.03 
18.32±0.04 
17.67±0.02 
17.50±0.01 
16.23±0.00 
18.73±0.03 
17.92±0.03 
17.12±0.01 
17.68±0.01 
18.()4±0.04 
16.76±0.01 
18.56±0.03 
17.14±0.01 
19.13±0.06 
16.!)6±0.01 
17.85±0.01 
16.08±0.00 
17.09±0.01 
18.37±0.02 
17.05±0.01 



15.70±0.01 
16.80±0.04 
15.44±0.01 
15.99±0.02 
16.02±0.01 
15.65±0.00 
16.06±0.01 
17.67±0.03 
17.84±0.04 
17.22±0.02 
16.99±0.04 
15.97±0.02 
18.40±0.16 
17.45±0.03 
16.72±0.01 
17.13±0.02 
18.53±0.06 
16.22±0.01 
18.24±0.05 
16.76±0.01 
18.46±0.06 
16.42±0.01 
17.25±0.02 
15.42±0.00 
16.64±0.01 
17.95±0.03 
16.41±0.01 



17.05±0.05 
17.79±0.08 
16.35±0.03 

16.57±0.03 
17.21±0.02 
16.38±0.01 
16.62±0.01 
17.71±0.05 
17.50±0.06 
17.33±0.04 
18.05±0.03 
17.47±0.01 
19.13±0.06 
18.29±0.07 
17.62±0.01 
18.26±0.02 
19.01±0.05 
17.23±0.01 
17.91±0.03 
17.19±0.01 
20.13±0.11 
17.09±0.02 
18.66±0.05 
16.90±0.02 
16.64±0.01 
18.06±0.05 
17.91±0.04 



15.59±0.02 
16.57±0.04 
14.89±0.01 
15.49±0.02 
15.84±0.01 
15.02±0.01 
15.42±0.01 
17.08±0.04 
16.62±0.04 
16.55±0.03 
16.85±0.02 
15.41±0.00 
18.08±0.04 
16.98±0.05 
16.82±0.01 
17.09±0.02 
18.43±0.04 
16.01±0.01 
17.12±0.03 
16.40±0.01 
19.04±0.07 
15.91±0.01 
17.47±0.02 
14.96±0.01 
15.81±0.01 
17.08±0.03 
16.47±0.01 



15.12±0.02 
15.97±0.06 
14.32±0.02 
15.07±0.03 
15.30±0.01 
14.39±0.01 
14.88±0.01 
16.61±0.05 
16.09±0.04 
16.07±0.03 
16.36±0.06 
15.92±0.01 
17.85±0.15 
16.58±0.04 
16.43±0.02 
16.56±0.02 
18.02±0.06 
15.57±0.01 
16.70±0.05 
16.09±0.02 
18.38±0.10 
15.38±0.01 
16.89±0.03 
14.33±0.01 
15.40±0.01 
16.64±0.05 
15.83±0.01 



0.051 
0.049 
0.050 

0.051 
0.034 
0.038 
0.049 
0.044 
0.047 
0.048 
0.046 
0.050 
0.044 
0.049 
0.043 
0.045 
0.041 
0.040 
0.038 
0.036 
0.037 
0.023 
0.025 
0.026 
0.025 
0.027 
0.025 



^The uncertainties included in the table are the intrinsic uncertainties associated with each measurment. 
Additionally, there is a systematic zeropoint uncertainty of 0.06 mag in the B-band, 0.05 mag in the Rc-band, 
and 0.05 mag in the Ic-band. 

^Photometric zeropoints for field Ul were determined relative to field N9. The zeropoint uncertainty for 
this field is 0.02 mag larger in each band than the other fields. See Section 2.2.2 for details. 



calibration of 0.3 mag. 

2.3. Spectroscopy 
2.3.1. Spectroscopic observations 

Spectroscopy was obtained using the Marcario Low Res- 
olution Spectrograph (LRS; Hill et al. 1998) on the 9.2m 
HET. Table 2 details the instrumental setup. The weather 
ranged from photometric to some cirrus. Typical image 
quality was ~2.6" FWHM during the 1999-2000 season, 
and ~2.3" FWHM during the 2000-2001 winter, measured 
on images taken before the spectrum was observed. Ad- 
ditionally, the shape of the point spread function (PSF) 
was much more regular and symmetric during the second 
season, as a result of image quality improvements made 
to the telescope. The only major instrumental change be- 
tween the two observing seasons was an upgrade of the 
CCD controller electronics from "version 1" in 1999-2000 
to "version 2" in 2000-2001. The version 2 controller pro- 
vides a much more stable bias and overscan level, faster 
readout, and lower readout noise. 

The typical exposure time was 30 minutes, and multiple 
observations were made for most of the galaxies to increase 
the signal-to- noise ratio (S/N) and aid in cosmic ray re- 
jection. Table 5 details the spectroscopic observations. 

In total we observed 19 LSB galaxies, two HSB ellip- 
tical galaxies (NGC 3872 and UGC 3844) and one HSB 
post-starburst irregular galaxy (NGC 1569). The HSB 
observations provide comparisons of old and young stellar 
populations, respectively. We also observed a K giant star 
(HD 12402) to use as cross-correlation template for the old 
stellar populations. 



2.3.2. Basic reductions of spectroscopic observations 

The spectral data reduction was done using IRAF. The 
observations spanned about 2.5 years beginning with early 
science operations of the HET, continuing through the 
telescope and instrument shake-down periods, and into 
normal operations. Consequently, we will go into some 
detail about the basic data reductions, and steps taken to 
find and correct systematic effects which might be present 
in the data. Readers not interested in the detailed steps of 
the basic reduction may skip ahead to Section 2.3.3. The 
steps of the basic data reduction are as follows: 

1) bias and overscan subtraction 

2) dark subtraction 

3) flatfielding 

4) cosmic ray cleaning 

5) wavelength calibration and rectification 

6) sky subtraction 

The dispersion axis in the LRS runs along rows. The 
version 1 CCD controller electronics, in use from 1999 Oc- 
tober - 2000 March, had a variable overscan level, depend- 
ing on the flux hitting the CCD. If there were rows with 
high flux such as the spectroscopic trace for bright objects, 
then the overscan level for those rows would be depressed. 
Consequently, for bright traces, the overscan had to be 
subtracted row by row. For faint traces, or uniform il- 
luminations such as with the flatfields, the overscan level 
was a smooth function of row number, and we subtracted 
a fit to the overscan, rather than doing it row by row. The 
version 2 controller produces a very stable overscan level, 
which is a constant for all rows regardless of the illumina- 
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Fig. 1. — Comparison between our photometry and that of OBC97 and OBSCI97. Both datasets have been corrected for galactic extinction. 
Sohd lines denote lines of equality; Dashed lines show the best fit to the zeropoint offset. Typical intrinsic uncertainties for all of the data are 
less than 0.05 mag, comparable to the size of the boxes, (a) Comparison of B-band magnitudes measured inside apertures of radius 5" . The 
fit to the zeropoint offset gives ^Tfizero — — 0.52 it 0.04 ina,g (rms — 0.19 mag). The two furthest outliers axe labelled and were not included 
in the fit. (b) Comparison of B-band total magnitudes. The zeropoint offset is l^vrizeLro = —0.11 ± 0.04 mag (rms=0.21 mag), excluding the 
two outliers N9-2 and Ul-4 (off the plot). The complex surface brightness profiles of these two galaxies were not well fit by OBC97 using 
single exponential disk fits, (c) Comparison of I-band 5" radius aperture magnitudes. Excluding N9-1 and N9-2, the mean zeropoint offset is 
Amzero = —0.31 ± 0.08 mag (rms=0.34 mag), (d) The I-band total magnitude comparison. Excluding N9-1, C5-3, and Ul-4 (off the plot) 
gives a best-fit zeropoint offset of Am^ero = 0.23 ± 0.08 mag (rms = 0.33 mag), (e) Comparison of (B— 1^) aperture colors. We derive a 
zeropoint offset of A(B— Ic)(5") = —0.24 ±0.07 mag (rms=0.29 mag). The two labelled outliers were not included in the fit. (f) Compajrison 
of the total (B— Ic) colors. The zeropoint offset is the same as for the aperture colors: A(B— Ic)(5") = —0.24 ± 0.07 mag (rms=0.30 mag), 
excluding the two labelled points. 



tion. 

Mean bias frames were produced each month by com- 
bining all bias (zero exposure time) images taken on nights 
when science data was obtained. These images were over- 
scan subtracted, then averaged together, with sigma clip- 
ping to eliminate cosmic rays. The coadded bias image 
was then fit with a low-order surface, and this fit was sub- 
tracted from all the science data for that month. The mean 
of the fit was about 1.4 ADU. Month to month variations 
were 0.1 ADU. 



Flatfield correction was done using lamp flats to cor- 
rect small scale pixel-to-pixel variations, and using twilight 
spectra to correct for the slit illumination in the spatial di- 
mension. Mean lamp flats were derived for each night, and 
checked against each other for variations, which were neg- 
ligible. While fringing is evident redward of ~7000 A, it 
is not seen to vary in the flats over the course of a month. 
There is some variation between months, caused by shifts 
in the slit position with respect to the image field center, 
and the corresponding wavelength range hitting the CCD. 
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Table 5 

IIET Spectroscopk' Ol!ser\:\ti()\s 



Name Exposure Time'^ Nobs^ PA slit*^ Dates of observations 



Pl-3 


1195 


1 


46 


Nov 30, 2000 


P6-1 


4200 


3 


45 


Dec 21,22(x2), 2000 


P9-4 


3600 


2 


41 


Nov 23,24, 2000 


Cl-2 


2014 


2 


-66 


Dec 10, 1999(x2) 


Cl-4 


5160 


3 


-67 


Nov 11, 1999(x2):Mar 11, 2000 


C3-2 


3300 


2 


-67 


Nov 10. 1999; Dec 7, 1999 


C4-1 


2770 


2 


-67, 67 


Nov 21, 2000: Jan 26, 2001 


C4-5 


3750 


2 


-67, 68 


Nov 12, 1999; Mar 6, 2000 




O ( uu 




-Do, Ol 


i\ov lyyy, uec ±4,1 (, ±yyy, iviar 10, zyjyjyj 


C5-5 


1666 


1 


-69 


Nov 30, 2000 


C6-1 


1200 


1 


-67 


Nov 24, 2000 


C8-3 


4800 


4 


-67, 67 


Dec 22, 2000(x2); Jan 27, 2001(x2) 


N3-1 


3300 


2 


-49, 47 


Dec 10, 1999; Mar 12, 2000 


N9-1 


3600 


2 


-75 


Nov 13, 1999; Dec 8, 1999 


N9-2 


5200 


3 


-76, 76 


Dec 13, 1999; Mar 3, 2000 


NlO-2 


3600 


2 


-66 


Dec 6, 1999; Dec 7, 1999 


NlO-4 


5280 


3 


-66, 66 


Dec 11, 1999; Mar 12,13, 2000 


Ul-4 


3600 


2 


-62 


Dec 13,14, 1999 


Ul-8 


1800 


2 


-60 


Dec 8,10, 1999 


NGC 3872 


900 


1 


-59 


Dec 7, 1999 


VGC 38-11 


Gon 


1 


-82 


Doc 7. 1909 


Ncx: i:)(ii) 


(iOU 


1 




No\- in, 200(1 



*^Total exposure time for all observations, in seconds 
'^Number of individual exposures 

"^Measured in degrees North through East; always parallactic 



The lamp flats taken for each month were combined into 
a single master flat. The rms uncertainty in these master 
flats is about 0.1% redward of 6000 A, increasing slowly 
blueward to 0.15% at 4800 A. 

Spectra of the sky taken during morning twilight were 
used to determine the slit illumination. The twilight flats 
were averaged together, then collapsed along the spectral 
dimension leaving a one dimensional "image" of the slit. 
This spatial profile was then applied to the lamp flat to cor- 
rect for the large scale response of the instrument, along 
the slit. The way the HET operates, the spherical pri- 
mary mirror is not moved during an observation. Instead, 
tracking is accomplished by moving the entire Prime Fo- 
cus Instrument Package (PFIP) relative to the primary 
mirror much in the way the Arecibo 305m radio telescope 
is pointed. As objects are tracked across the mirror, the 
pupil entering the instrument changes and, consequently, 
so does the slit illumination. We found that the twilight 
flats taken with the PFIP location fixed often did not rep- 
resent the slit illumination during our science exposures. 
The gross vignetting pattern of the telescope optics was 
represented, but small gradients remained. These gradi- 
ents were corrected as part of the sky subtraction proce- 
dure, described below. 

Cosmic ray rejection was done using an IRAF script 
called speccrrej (a prototype of a task which is planned 
for the Gemini IRAF package). This task makes a two 
dimensional model of the longslit spectrum, sensitive to 
features larger than ~ 30 square pixels. It does a fairly 
good job of simulating the spectrum, as well as most of 
the sky lines. However, sharp features are not modeled, 
and show up in the difference image. Most of these fea- 
tures are cosmic rays, though strong skylines and some- 
times strong emission lines or deep absorption troughs are 
also flagged. The cosmic ray masks were all checked man- 
ually to make sure that real spectral features were not 
flagged. The flagged cosmic rays were then replaced by 



the values from the model spectrum. Pixels affected by 
cosmic rays were added to the bad pixel masks for each 
exposure. 

Spectral arc lamps (Ne and HgCdZn) were observed 
each night with our setup. We used the IRAF tasks 
identify, reidentify, and f itcoords to establish the 
spectral dispersion function in the two dimensional spec- 
tra. There were 34 lamp lines used over the full wavelength 
range, and the fifth order polynomial fits typically had an 
RMS scatter of 0.2 A. There were small (less than 2 A) 
wavelength offsets found for the skylines in many of the 
science exposures, caused by a shifting of the slit position 
between setups. These shifts caused a wavelength shift in 
the dispersion function, but no change in its shape. The 
wavelengths of several strong skylines (OIA5577, OIA6300) 
were used to correct for these offsets. The identify task 
was also used to find the spectral trace along the spatial 
dimension. In cases where the trace of a galaxy was too 
faint to be seen, the fit for a brighter object observed the 
same night was used instead. We then rectified the two 
dimensional spectra. 

Prior to sky subtraction, we checked the quality of the 
spatial illumination correction for each spectrum, by mea- 
suring the fiux of several strong sky lines along the slit. 
The sky level should be constant along the slit, but due 
to the varying illumination of the HET image pupil, the 
simple twilight fiatfield was not always effective. In these 
cases, we fitted any residual illumination gradient with a 
smooth low order function, and divided it out. Sky regions 
bracketing the spectral trace were averaged, and then sub- 
tracted from the full 2D spectrum. This produced good 
results, except near the OIA5577, NaDA5895, and OIA6300 
skylines, which were simply interpolated over. These bad 
regions are marked in the spectra presented in Section 3. 

2.3.3. Spectral extraction, coaddition, flux calibration 
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We extracted the spectra in apertures of length rg (thus 
including the signal within re/2 of the center, rounding 
off to the nearest integer pixel), by simply summing the 
appropriate rows. Since the galaxies cover a wide range 
of distances, the fixed slit width of 2" will cover a vary- 
ing fraction of the galaxy disks. However, since no good 
spectral gradient data exist for these types of galaxies, we 
choose to simply note this fact and not attempt any cor- 
rections for it. The galaxies C3-2 and C4-1 did not have 
exponential scale lengths reported in OBC97. Wc chose 
the apertures for these galaxies based on their light pro- 
files along the slit and to optimize the S/N; we use an 
aperture length of 6" for C3-2, and 11" for C4-1. 

A special correction had to be made for the data from 
1999 November. During the first month of science ob- 
servations with the telescope, a glass membrane pellicle 
was used as part of the telescope guiding system. This 
membrane caused a fringe pattern to appear in the data, 
with a fringe wavelength of 300-400 pixels, and a strength 
and location that varied from one exposure to the next. 
For each of the galaxies observed in 1999 November, we 
obtained at least one additional spectrum during a later 
month. To remove the fringe pattern from the data, we 
fit the continuum, and normalized the spectrum by this 
fit. We then multiplied the normalized spectrum with the 
average continuum fit from the data not affected by the 
fringing problem. 

Most of the galaxies had multiple observations, often 
taken several months apart. Before coadding the spectra, 
we corrected the wavelength scales to a consistent velocity 
with respect to the sun. The spectra were then averaged. 

Each night of science observations, a spcctrophotomct- 
ric standard star was observed. The relative flux response 
of the instrument is very stable (with the exception of the 
data from 1999 November), and we based the relative flux 
calibration on all the standard stars observed during an 
observing season. For the 1999-2000 season, wc used 15 
standard star observations, and for the 2000-2001 season 
we averaged 7 observations. The stars observed were: HD 
84937, HD 19445, HZ 4, HZ 44, Feige 34, Feige 66, Feige 
67, Gl 191B2B, Hiltner 600, BD+25 3941, and BD+26 
2606. The observed spectral energy distributions (SED) 
were scaled to account for grey extinction due to cirrus, 
and then fit with a sixth order polynomial, with an RMS 
deviation of less than 0.1%. When the different observa- 
tions are scaled to the same value at 5800 A, the full range 
of measurements at both 4800 A and 7150 A spans 0.25%, 
which we take to be the maximum expected relative un- 
certainty over large wavelength ranges. We then applied 
these relative flux calibrations to all the coadded spectra. 

In parallel with the spectral reductions and extractions, 
variance spectra were produced. These are derived by first 
assuming poisson noise from the sky and object counts, 
then folding in readout noise, flat-fielding uncertainty, the 
effects of rectification and extraction, and finally the co- 
addition of spectra. 

2.3.4. Redshifts 

Redshifts are derived in two separate ways, (1) from 
cross-correlation with absorption line templates and (2) 
from measurements of the emission line wavelengths. 

We used the IRAF task xcor to compute the cross- 
correlations for galaxies showing stellar absorption fea- 



tures, with an observation of HD12402 (spectral type 
Kini) as the stellar template. The restframe wavelength 

range used for the computation was 4300 A - 6500 A 
for galaxies with no significant H/3 emission (Cl-2, Cl- 
4, Ul-8, N3-1). For galaxies with H/3 emission, the cross- 
correlation was performed over the range 5000 A - 6000 A, 
which still includes several strong absorption features. The 
uncertainties determined by the xcor task were 60 km/s 
for the galaxies without H/3 emission and 75 km/s for the 
systems with both absorption and emission features. 

Redshifts for the emission line galaxies were based 
on positions of the lines H/3, [O III]AA4959,5007, Ha, 
[N II]A6584, and [S II]AA6717,6731. We take the mean ve- 
locity of all the strong lines as the systemic velocity, and 
the rms scatter as the uncertainty in the velocity. This 
uncertainty was typically ^ 30km/s for galaxies with at 
least four emission lines detected. 

For the galaxies where we measiirc redshifts from both 
the emission lines and stellar absorption features, the two 
measurements were always consistent with no velocity off- 
set, to within the uncertainties. We transform the red- 
shifts into the heliocentric reference frame. The measured 
redshifts are listed in Table 7. We list the emission line 
redshifts when available, otherwise the absorption line red- 
shifts are shown. 

2.3.5. Emission line measurements 

Emission line strengths are measured using the IRAF 
task f itprof s. Nearby lines are measured simultaneously, 
using the deblend option, and constraining the fit to have 
a single width for all the nearby lines. We use the variance 

spcKjtra to estimate the "inverse gain" parameter, used by 
f itprof s in its Monte Carlo estimation of the uncertain- 
ties. 

In many cases, the measurement of the H/3 emission 
strength is clearly affected by underlying stellar absorp- 
tion. We correct this individually for each galaxy. The 
adopted correction procedure is as follows. The wave- 
length coverage of the emission line is defined at its base, 
the widest point. The continuum level is defined with 
a straight line between the two sides of the base. The 
emission line flux and equivalent width are then measured 
above this continuum level. The V1999 single burst stel- 
lar population model spectra, convolved to our instrumen- 
tal resolution, are then used to estimate the effect of the 
underlying stellar absorption profile on the emission mea- 
surements. For each galaxy, we find the model which best 
fits the absorption lines in the restframe wavelength range 
4800A < A < 5450A. We also determine what range of 
model ages and metallicities brackets a reasonable fit to 
the spectrum. The fits are all done by eye. We define 
a continuum level on the models, by drawing a straight 
line between the model values at either side of the emis- 
sion line wavelength interval, and measure the equivalent 
width of absorption below this level. We then combine 
the equivalent width of absorption seen in the models with 
the equivalent width of emission measured above the in- 
terpolated continuum in the data to get a corrected H/3 
equivalent width measurement. Finally, we scale the H/3 
flux measurement by the ratio of the corrected to uncor- 
rected equivalent widths to determine the total emission 
line flux corrected for underlying absorption. Examples 
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Fig. 2. Correction of the H/3 emission line strengths for underlying H/3 absorption. Solid lines - observed spectra; dotted lines - best fitting 
model; dashed lines — youngest acceptable model; dot-dashed lines - oldest acceptable model. The strongest emission and absorption features 
in this wavelength region are labelled in panel (a), (a) C6-1 with the best fitting model, agc=3.2 Gyr, [Fe/H]=-0.4. (b) C6-1 with oldest 
(age=13.2 Gyr, [Fe/H]=-0.7) and youngest ( age=1.0 Gyr, [Fe/H]=0.0) acceptable models, (c) C4-1 with the best fitting model, age=250 
Myr, [Fe/H]=0.0. (d) C4-1 with oldest (age=630 Myr, [Fe/H]=0.0) and youngest (age=100 Myr, [Fe/H]=0.0) acceptable models. Note that 
the depth of the H/3 absorption reaches a maximum for an SSP of 200 Myr, and so appears weaker in both the older and younger models. 



of this technique appHed to two of the galaxies are shown 
in Figure 2. Table 6 lists the measured uncorrected H/3 
equivalent widths as well as the parameters of the best fit- 
ting, youngest, and oldest models, and their corresponding 
absorption equivalent width corrections. 

We list two sets of uncertainties for the H/3 fluxes in Ta- 
ble 7. The first set are the random uncertainties derived 
from the fit to the emission line flux, the same as for the 
other emission lines. Additionally, we include a set of un- 
certainties showing the range of the absorption corrections 
consistent with the reasonable fitting models. 

The V1999 models, which use an empirical spectral li- 
brary, are limited in their ability to model young ages 



at sub-solar metallicities. We also investigated using the 
StarburstQQ models (Leitherer et al. 1999), which cover 
very young ages at both solar and 1/20 solar metallicity. 
However, the narrow wavelength coverage of their models 
does not permit us to fit the metallic absorption features 
at longer wavelength, which provide the most constraint 
on the range of reasonable models. 

Telluric absorption in the atmospheric B band (A ~ 6840 
A) affected our line strength measurements for two galax- 
ies. The [N II]A6584 line was absorbed in the galaxy C5-3. 
Nebular physics predicts [N II]A6584 = 2.95 ■ [N II]A6548 
(Osterbrock 1989, p. 61). For the 9 LSB galaxies with 
emission and where strong Ha does not make it impossi- 
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Table 6 

h/3 equivalent width corrections 









hnst-fit moHpl'^ 








min RWCH/1~I rr,ndr]'= 


Name 


EWCHB),.^ 


age 


[Fe/Hl 


EW(Hfl')„h, 


age 


[Fe/Hl 


EW(H/3'),h« 


age 


fFo/Hl 






(A) 


(Gyr) 




(A) 


(Gyr) 




(A) 


(Gyr) 




(A) 


C3-2 


2.0±0.1 


1.58 


0.2 


-1.3 


0.63 


0.0 


-1.0 


11.0 


-0.7 


-1.9 


C4-1 


7.5±0.6 


0.25 


0.0 


-2.2 


0.10 


0.0 


-1.8 


0.63 


0.0 


-2.2 


C4-5 


0.9±0.3 


0.79 


0.0 


-1.2 


0.40 


0.0 


-1.0 


1.58 


0.0 


-1.5 


06- 1 


2.7±0.3 


3.16 


-0.4 


-1.6 


1.00 


0.0 


-1.3 


13.2 


-0.7 


-2.4 


C8-3 


3.4±0.4 


0.50 


0.0 


-1.1 


0.10 


0.0 


-0.9 


0.79 


0.0 


-1.2 


NlO-2 


1.3±0.1 


7.94 


-0.4 


-1.0 


1.00 


0.2 


-0.8 


17.4 


-0.7 


-1.4 


N9-1 


11.5±0.3 


0.40 


0.0 


-2.6 


0.03 


0.0 


-1.7 


1.00 


0.2 


-2.9 


N9-2 


27.7±0.9 


0.25 


0.0 


-4.4 


0.10 


0.0 


-3.4 


0.63 


0.0 


-4.5 


Pl-3 


7.2±0.4 


1.00 


0.0 


-1.3 


0.40 


0.0 


-1.0 


6.31 


-0.7 


-1.9 


P6-1 


1.3±0.2 


3.16 


-0.4 


-1.0 


1.00 


0.2 


-0.8 


11.0 


-0.7 


-1.4 


Ul-4 


2.6±0.1 


0.63 


0.0 


-2.0 


0.50 


0.0 


-1.9 


0.79 


0.0 


-2.3 



^Parameters for the "best-fit" stellar population model and the equivalent width of H/3 absorption over the narrow wavelength region 
for which we measure emission. The "best lit" as well as the determination of the maximum and minimum reasonable corrections 
were done by eye. See Section 2.3.5 for details. 

''The model which requires the largest correction to the H/3 equivalent width, but is still a reasonable fit to the other absorption 
features in the spectrum. This is typically the youngest, highest metallicity, marginally acceptable model. 

'^The model which requires the smallest correction to the H/3 equivalent width, but is still a reasonable fit to the other absorption 
features in the spectrum. This is typically the oldest, lowest metallicity, marginally acceptable model. 



ble to measure the [N II]A6548 line, the mean [N II]A6584 / 
[N II]A6548 ratio is 2.9±0.9. Wc use the measured value of 
[N II]A6548 in C5-3 and the theoretieal ratio to ealeulate 
the non-absorbed [N II]A6584 strength (0.27±0.08 relative 
to H/3). For the galaxy N9-1, [N II]AA6548,6584, and Ha 
all suffered some absorption. For this galaxy, we estimate 
the additional uncertainty on the [N II]A6584/Hq ratio to 
be 10%, and we also likely underestimate the Ha/H/3 ra- 
tio by 10%. The uncertainties on the relative emission line 
fluxes listed for N9-1 in Table 7 have been adjusted to ae- 
count for this. In the following the [N II]A6584/Hq; ratio 
is simply called [N II] /Ha. 

We correct the relative emission line strengths for galac- 
tic reddening using the Cardelli extinction law (eqs. 3a and 
3b in Cardelli et al. (1989)), with the galactic reddening 
measurements from Schlegel et al. (1998) as listed in Table 
1. Ry=3.1 is assumed. The amount of galactic reddening 
is variable across the field of NGC 1569. We follow Kob- 
ulnicky & Skillman (1997) and adopt the value of 0.5 mag 
(Burstein & Heiles 1982), so that we may make a direct 
comparison to their results (see section 2.3.7). We have 
not corrected the emission line strengths for dust redden- 
ing internal to the galaxies. We discuss the dust content of 
these galaxies in Section 7. The emission line ratios used 
for ionization determination and abundance analysis (e.g. 
[O III]A5007/H/3 and [N II]/Ha) utilize lines with small 
wavelength separations, and these ratios are not affected 
by reddening at the level seen in these galaxies. The Ha 
equivalent widths are also corrected for the effects of cos- 
mological expansion, EW(Arest)=EW(Aobs) ' + z) (see 
Peterson 1997, page 156). 

We present the measured emission line strengths as well 
as the equivalent width of Ha in Table 7. These data 
do not have an absolute flux calibration. There are sev- 
eral reasons for this. As objects are tracked across the 
sky by moving the HET prime focus instrument platform, 
the region of the primary mirror seen by the instrument 
changes. Furthermore, not all of the primary mirror is 
seen all the time, and thus the telescope throughput varies 
during an object track. Additionally, the PSF at the spec- 



trograph entrance slit changed (sometimes dramatically 
increasing in size) during exposures, and some of the obser- 
vations were taken through cirrus clouds. Consequently, 
we tabulate only relative fluxes, and for convenience we 
have scaled them to an H/3 flux of 1. For Table 7, we 
list both the random and systematic uncertainties on the 
scaled H/3 flux, rather than incorporating these into the 
uncertainties of the other flux measurements and listing 
line ratios. The uncertainties on H/3 are incorporated into 
the [O III]A5007/H/3 ratios plotted in Figure 6 and the 
Ha/H/3 ratios shown in Figure 11. 

2.3.6. Absorption line indices 

We have measured the strength of the stellar absorption 
lines H/3, Mgb, Fe5270, and Fe5335 using the line index 
definitions of the Lick/IDS system (Worthey et al. 1994). 
Several calibrations were made to match our observations 
to the characteristics of the IDS spectrograph and put our 
measurements onto the Lick system. Before measuring 
the indices, it was necessary to match our spectra to the 
instrumental resolution of the IDS spectrograph. We con- 
volved our spectra with a wavelength dependant Gaussian 
to accomplish this. Additionally, stellar features in galaxy 
spectra are broadened by the internal velocity dispersion 
of the galaxies, and this broadening will have an effect on 
the measured strength of the line indices. To facilitate 
comparison with stellar population models, we correct our 
observations to the equivalent strength for a system with 
zero velocity dispersion, following the method of Davies, 
Sadler, & Peletier (1993): we determined the correction 
by taking a spectrum of a K giant star, already convolved 
to the IDS resolution, and convolving it with Gaussians 
of increasing velocity width, in 25 km/s increments, up to 
400 km/s. We measured the index strengths on this suite 
of stellar spectra, and used the offsets between each ve- 
locity width and the unbroadened spectrum to correct the 
galaxy measurements. Since we do not actually resolve 
the velocity dispersion of the LSB galaxies, we assume a 
value of 75 km/s. The correction factors by which the 
measurements are multiplied are as follows: H/3 (1.004); 



Table 7 
Spectroscopic Data 



to 



Name 






M 


c 

B 


H/3 


Mgb 




< Fe > 


EW(Hcv ) 






[OIII] A5007 


A6548 


Ha 


[NII]A6584 


[SII]A6717 


[SII] A6731 


Pl-3 


19 


17505 


-19 


.67 




0.9±0, 


.4 


1.4±0 


.4 


31±1 


1 


.nn±o 


„- + 0.07 
■ —0.04 


1.48 + 0.05 


0.20 + 0.04 


3.44+0.07 


0.70+0, 


06 


0.89 + 0.05 


0.65 + 0.05 


P6-1 


34 


10891 


-18 


,53 




3.1±0 


.2 


2.1±0, 


.2 


8.0±0.2 


1 


,00±0 


l^^O 




0.75 + 0.07 


5.0 + 0.1 


2,3 + 


1 


1.70+0.09 


1.13 + 0,08 


P9-4 
Cl-4 


7 

70 


4108 
3883 


-14 
-17 


,67 
,35 


1.9±0.1 


3,5±0 


1 


2,2±0 
2,50±0 


3 

09 


73±7 


1 


.00±0 


06 


3.1 + 0.1 


0.04 + 0.04 


3.76 + 0.0.09 


0.30 + 0. 


05 


0.56 + 0.06 


0.36 + 0.06 


C3-2 


41 


4528 


-17 


,63 




2,4±0 


2 


2,0±0 


2 


ll.OiO.l 


1 


,00±0 


+ 0,18 
""-0,09 




0.33 + 0.04 


3.84+0.05 


1.49+0, 


04 


0.82 + 0.05 


0.76 + 0.06 


C4-1 


14 


4927 


-16 


,48 




1,0±0 


5 


0,8±0 


5 


30±2 


1 


,00±0, 


07+"-™ 
— 0,04 


1.53 + 0.07 


0.09 + 0.04 


2.64 + 0,09 


0.25 + 0, 


05 


0.56 + 0.06 


0.39 + 0.05 


C4-5 


23 


4633 


-16 


,01 




1.1±0 


.3 


1.7±0 


.3 


6.3±0.3 


1 


,oo±o. 


2^ + 0,12 
^"-0,13 


0.5 + 0.1 


0.47+0.13 


2.9+0.1 


0.9 + 


.1 


1.3+0.1 


0.9+0.1 


C5-3 
C5-5 


12 
5 


12944 
5473 


-17 
-16 


.21 
.49 




0.9±0 


.6 


1.0±0 


.5 


63±1 
25±2 


1 
1 


.00±0. 
.00±0. 


06 
28 


3.78 + 0.09 
4.3 + 0.4 


0.09+0.03 


3.53+0.08 


0.27+0. 
5.8+0 


08*^ 
.4 


0.62+0.05 


0.41+0.04 


C6-1 


22 


25035 


-20 


.13 




3.1±0 


3 


2.2±0 


3 


14.4±0.6 


1 


.00±0. 


12+0.17 
-0.08 




0.65+0.06 


4.3+0.1 


1.96+0. 


09 






C8-3 


17 


3648 


-15 


.91 




1.1±0 


4 


1.0±0 


4 


16.3±0.6 


1 


.00±0. 


,,+0.03 
^"^-0.03 


3.2+0.2 


0.25+0.06 


3.2+0.1 


0.31+0. 


07 


0.74+0.09 


0.42+0.07 


N3-1 


72 


25222 


-19 


.67 


1.7±0.1 


4.7±0 


1 


3.03±0 


09 






















N9-1 


38 


15128 


-19 


.38 




1.6±0 


2 






41±1 


1 


.on±o. 


02 + 0.02 
"^-0.07 


0.71 + 0.02 


0.12 + 0.01 


2.88+0.02 


0.80+0. 


01 


0.81 + 0.02 


0.57+0.02 


N9-2 


13 


73629 


-19 


.84 




1.3±0. 


6 


1.3±0. 


5 




1 


.00±0. 


n-i+o.oo 

"^-0.04 


1.09 + 0.03 














NlO-2 
NlO-4 


38 
8 


20668 
7447 


-20 
-16 


.69 
.20 




2.8±0. 


.2 


2.2±0. 


.2 


8.1±0.2 
7±1 


1 


.00±0. 


11+0.17 
-0.07 


0.33 + 0.08 


0.70 + 0.09 


3.80+0.08 


1.90+0. 


06 


"■^'^+0.07 


32t°-^' 
"•"^±0.07 


Ul-4 


53 


3416 


-16 


.86 




1.0±0. 


.2 


0.8±0. 


.1 


11.9±0.2 


1 


.00±0. 


"^-0.03 


1.13 + 0.04 


0.11 + 0.02 


2.18+0.03 


0.32+0. 


03 


0.57+0.02 


0.43+0.02 


Ul-8 


35 


32123 


-20 


.27 


2.0±0.3 


3.0±0. 


.2 


2.2±0. 


.2 






















NGC3872 
UGC3844 


153 
38 


3114 
3192 


-21 
-20 


.838 
.778 


1.47±0.05 
1.5±0.2 


5.04±0. 
4.6±0. 


.05 
.2 


2.83±0. 
3.5±0. 


.04 
.2 
























NGC1569'' 


120 


-96 


-17 


.238 




0.68±0. 


,07 


0.71±0. 


.06 


95.6±0.6 


1 


.00±0. 


00 


5.96 + 0.01 


0.06 + 0.00 


3.92+0.00 


0.14+0. 


00 


0.20+0.00 


0.14+0.00 



^Signal- to-noise per A of the continuum, measured between 5050A and 5150A. 

^Heliocentric radial velocity: Baaed on emission lines, if present. Otherwise based on cross-correlation with absorption line template. 

*^Total absolute magnitude: Based on our photometry, if available, otherwise from OBC97, Assuming Ho=75 kms ~ Mpe ~ , and distances based on redshifts in the CMB reference frame. 

•^Emission line fluxes are normalized to f(H^)^l. The first set of uncertainties on H/3 are the random uncertainties and the second set are the range of absorption corrections consistent with the 
stellar population models. If no second set is listed, the correction is negligible compared to the uncertainty due to random error. 

^Calculated from [NII]A6548 X 2.95 

^Lower value is the intrinsic uncertainty, upper value is a systematic uncertainty, caused by poor flux calibration at the extreme red end of the spectrum. See Figure 6. 
^Data from the LEDA database (http://leda.univ-lyonl.fr), and adjusted to Ho^75 kms~ Mpc~ . 
^Uncertainties listed as 0.00 are less than 0.005. 
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Mgb (1.02); Fe5270 (1.02); Fe5335 (1.04). The correction 
is such that for galaxies whose true velocity dispersion is 
larger than we have assumed (75 km/s), we will derive 
index values which are weaker than the actual values for 
that stellar population if it were observed with zero intrin- 
sic velocity dispersion. If the actual velocity dispersion of 
the galaxy is 150 km/s instead of 75 km/s, then we will 
systematically underestimate the H/3 strength by 1%, the 
Mgb strength by 3%, the Fe5270 strength by 5%, and the 
Fc5335 strength by 8%. For N3-1 and the two elliptical 
galaxies we correct for a velocity dispersion of 250 km/s. 
NGC 3872 has a velocity dispersion of f7=243km/s (Davies 
et al. 1987), while we simply assume 250 km/s for UGC 
3844 and N3-1. The correction factors arc H/3 (1.02); Mgb 
(1.17); Fc5270 (1.16); Fe5335 (1.36). 

Uncertainties for the index strengths are derived from 
the variance spectra. Essentially, we use poisson statistics 
on the flux in both the line and continuum bandpasscs, 
additionally accounting for the contribution caused by the 
readout noise, and flatfielding, relative flux calibration, 
and sky subtraction uncertainties. 

In several cases, residuals from imperfect subtraction 
of strong skylines (OA5577, NaD) fall in the index band- 
passes. When the residuals lie in the index continuum 
bandpasses, we correct the residuals by interpolating over 
the affected region, using values determined from compar- 
ison with the V1999 stellar population models. Observa- 
tions corrected in this way arc NlO-2 (Mgb red contin- 
uum, Fe5270 blue continuum), N3-1 (Mgb blue contin- 
uum), Ul-8 (Fe5270 red continuum), and C6-1 (Mgb blue 
contimmm). In some cases, the sky residuals lie on the 
absorption line itself, and this cannot be corrected. Both 
Fe5270 and Fe5335 arc affected in N9-1. The Fe5335 line 
is affected for C5-3 and Ul-8 while the Fe5270 hne is af- 
fected for Pl-3. There are strong correlations between 
Fe5270, Fe5335, and <Fc> (defined as <Fe> = (Fe5270 
-|- Fe5335) / 2). Using the LSB galaxy data not affected 
by the large sky subtraction uncertainties, we derived the 
relation Fe5335 = (1.14±0.08)- Fe5270 -(0.74±0.24). We 
use this relation to correct the affected indices. 

Intrinsic emission from these galaxies could affect the 
line indices Mgb and H/?. Goudfrooij & Emsellem (1996) 
have found that the [NI] emission doublet at 5199 A can 
affect the strength of the Mgb index in early-type galax- 
ies which have ionized gas emission. Many of these LSB 
galaxies have emission lines, so this is cause for concern. 
However, the [NI] doublet is strong in LINERs, not H II 
regions, and none of these galaxies show line ratios consis- 
tent with this type of nuclear activity (see Section 4) . Ad- 
ditionally, since these are integrated spectra, rather than 
nuclear spectra, we do not expect the Mgb index to be 
contaminated. The H/3 index, on the other hand, is cer- 
tainly affected by emission in all but three of the galaxies 
in our sample. For one more galaxy (Ul-8), the 11/3 emis- 
sion must be small. The observed wavelength range for 
Ul-8 does not include either Ha or [O II]A3727. While 
the spectrum is dominated by absorption features, there 
is a small amount of emission seen at [O III]A5007 which 
suggests that the H/3 absorption index will also be affected 
by emission. The equivalent width of [O III]A5007 emis- 
sion is 1.2±0.4 A. Trager et al. (1999) determine the re- 
lation between [O III]A5007 and H/3 emission strength in 



eUiptical galaxies to be EW(H/3) ~ 0.6-EW([O III]A5007). 
This ratio is metallicity sensitive, and will tend to decrease 
with decreasing metallicity. The integrated galaxy spec- 
tra studied by Kobulnicky, Kennicutt, & Pizagno (1999) 
show ratios as low as EW(H/3) 0.25-EW([O III]A5007) at 
the metallicity appropriate for Ul-8 (~ 0.4 Z©; see Figure 
7). We therefore apply a correction to the H/3 absorption 
index for Ul-8 of 0.5 ± 0.3 A. 

For the 15 LSB galaxies with strong emission lines, a 
correction of the H/3 absorption index for the H/3 emission 
(and its uncertainty) would be a significant fraction of the 
H/3 absorption strength. Consequently, we do not present 
this absorption index for any of the other LSB galaxies. 

The absorption line strengths of H/3 , Mgb, and < Fe > 
are presented in Table 7. 

2.3.7. Comparison to literature data 

Our sample has a very wide range in radial velocities, 
from 3000 km/s to more than 60,000 km/s (z=0.2). These 
same galaxies have been observed in the H I survey of 
O'Neil et al. (2000a), who give detections for twelve of 
them. We find several discrepancies between our redshift 
measurements and the redshifts reported in that survey. 
The H I survey only covered the velocity space out to 
13,000 km/s so more distant objects would not have been 
detected. However, for the galaxies Pl-3, C6-1, and N9-2, 
lower H I redshifts must have been attributed incorrectly 
to these higher redshift galaxies. The H I detections for 
Pl-3 and C6-1 can be explained by beam confusion. At 
21cm, the beam-size of the Arecibo 305m radio telescope 
is almost 3', and both of these galaxies appear in the sky 
at a projected distance of less than 3'from a large spiral 
galaxy at lower redshift. The H I detection ascribed to Pl- 

3 is probably NGC 7631 (czhi = 3746 km/s; czopt = 3741 
km/s (Falco et al. 1999; hereafter UZC), and the detection 
ascribed to C6-1 is probably NGC 2595 (czhi = 4322 km/s; 
czopt = 4320 km/s, UZC). The mis-identification also ex- 
plains why these galaxies appeared to have extreme Mhi/L 
ratios, and lie more than 5a from the I band TuUy-Fisher 
relation (O'Neil et al. 2000a). Similar cases of beam con- 
fusion explain discrepant measurements for the lower red- 
shift LSB galaxies as well. The detection ascribed to C8-3 
is probably UGC 4308 (czhi = 3570 km/s; czopt = 3566 
km/s, UZC), and the detection of C5-5 is more likely UGC 
4416 (czhi = 5524 km/s; czopt = 5555 ± 45 km/s, UZC). 

There are also cases where the optical and H I redshifts 
are discrepant, but beam confusion is not obviously the 
reason: N9-2 (listed above), C4-1, Cl-2, and possibly P9- 

4 (where the optical and H I velocities differ by 3a). Cl-2 
in particular is a hard case to explain the H I detection. 
The difference between optical and H I velocities is ~3900 
km/s so the H I detection is clearly a different object. 
However, there are no obvious candidates within 15' of 
the pointing, which lie at the redshift of the H I detection, 
and very few bright objects in that area for which opti- 
cal redshifts were not obtained by the redshift survey of 
Zabludoff & Mulchaey (2000). We are forced to conclude 
that either the H I observation was errant, due to either 
noise or incorrect pointing, or that there are very large H I 
clouds with little associated optical light. Two starless ob- 
jects have been found by the HIPASS survey (Kilborn et 
al. 2000; Ryder et al. 2001) but the first one had only a 4 
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km/s velocity width, and the second was associated with 
the galaxy NGC 2442, so neither one is a good analogue. 

The remaining four galaxies, C5-3, P6-1, NlO-4, and 
Ul-4 are in good agreement, less than la separating the 
optical and H I measurements. 

Recently, Chung et al. (2002) used the VLA to map the 
H I around the galaxies Pl-2, Pl-3, C6-1, C4-1, and C4-2 
from the OBC97 sample. They also discovered the beam 
contamination issues affecting the O'Neil et al. (2000a) 
Arccibo 305m singic-dish results. The maps for Pl-2, Pl- 
3 and C6-1 show clearly that the H I is all at the position 
of neighboring spiral galaxies. The H I detection for C4-2 
was confirmed, but the map of C4-1 did not detect any H I 
even though the VLA observations were sufficiently deep 
to detect the H I, if it were present at the level measured 
at Arecibo. This result suggests that the N9-2 and Cl-2 
H I detections may also be noise artifacts, and not starless 
gas clouds. 

We compare our emission line strength measurements 
for the post-starburst dwarf galaxy NGC 1569 with the 
results of Kobulnicky & Skillman (1997). Our aperture 
lies perpendicular to the C and D slit placements of their 
work and near the two super-star clusters. Kobulnicky 
& Skillman found that the combination of galactic and 
intrinsic reddening is both large and variable over the 
face of this galaxy. In accordance with their methods, 
we choose a value for the total extinction (galactic plus 
intrinsic) which brings our measured Hq/H/3 ratio to 
2.83 after applying the reddening correction. We mea- 
sure a total E(B— V) of 0.83, which is within the range 
of values found by Kobulnicky & Skillman (following 
Kobulnicky & Skillman and Burstein & Heiles (1982), 
we assume that 0.5 mag of the reddening originates in 
our galaxy). Once we apply the extinction correction to 
all the emission lines, our relative flux measurements for 
the lines H7 (0.45±0.01), [O 111]AA4959,5007 (1.94±0.01, 
5.72±0.01, respectively), [S III]A6312 (0.008±0.001), 
[N 11]AA6548,6584 (0.042±0.001, 0.098±0.001, respec- 
tively), HeIA6678 (0.023±0.001), and [S II]AA6717,6731 
(0.137±0.002, O.lOOiO.OOl, respectively) all lie within the 
range covered by the Kobulnicky & Skillman measure- 
ments for this part of the galaxy. The only discrepant line 
is [O III]A4363 (0.124±0.007), which lies on the blue end 
of our spectrum where the flux calibration is less secure. 
Over the wavelength range 4800 A to 7000 A, our emission 
line measurements are fully consistent with the results of 
Kobulnicky &: Skillman. 

The absorption line index strengths for NGC 3872 are 
compared with the results of Trager et al. (1998). We 
extract the spectra using a spatial aperture equivalent to 
the Lick/IDS aperture, and make velocity dispersion cor- 
rections appropriate for (TngC3872 = 250 km/s. Trager et 
al. (1998) measure the following values: H/3: 1.28±0.22; 
Mgb: 4.91±0.30; <Fc>: 2.22±0.22. Our measurements 
for H/3 and Mgb are within la of their values. For <Fe>, 
our measurement is about 2.5a higher than theirs. 

3. PRESENTATION OF THE SPECTRA 

We present the spectra in Figures 3a-e and 4. The 
spectra for these 19 LSB galaxies qualitatively resemble 
the spectra of HSB galaxies over the full range of spec- 
tra seen for all Hubble types, from ellipticals all the way 
through late- type spirals and starbursting irregular galax- 



ies. For comparison, we show our spectra of the two ellip- 
tical galaxies, NGC 3872 and UGC 3844, in Figure 3a, and 
the spectrum of the post-starburst irregular galaxy NGC 
1569 in Figure 3d. 

The emission line galaxies have been ordered by the 
strength of their [N II]A6584 emission relative to Ha, 
from strongest to weakest. This is essentially a sequence 
of decreasing metallicity as we will show in Section 6. 
The strength of the [O III]A5007 emission line increases 
along the sequence, due to a combination of the decreasing 
metallicity and the associated decreasing amount of dust 
reddening in the emission regions. This will be discussed 
in Section 7. 

The spectral extraction was done using global apertures, 
which blend together stellar light and line emission from 
both H II regions and diffuse emission regions. While we 
do not generally have sufficient spatial resolution to ob- 
serve gradients in the galaxy properties, we note that the 
galaxy Ul-4 has a clear inner and outer region. The inner 
region of radius 3" is dominated by a bright continuum 
with stellar absorption features. Outside this region there 
is a peak in the emission on both sides of the nucleus at 
about 5". Further out, both the emission and continuum 
fall off. This is illustrated in Figure 5. The spectroscopic 
aperture we use for this galaxy is 2" x 32" and includes 
light from both the inner and outer regions. 

Analysis of the galaxies in Figure 3(a-d) will be pre- 
sented in the following sections. Here we discuss in detail 
the galaxies presented in Figures 3e and 4. 

The galaxy N9-2 (Figure 4) has a redshift of 0.246 
(cz=73629 km s"-'^). The restframe wavelength coverage 
for this galaxy is significantly bluer than the rest of the 
sample, and we do not have observations of the wavelength 
region near Ha. After applying the (1+z)"^ cosmological 
surface brightness correction and a fee correction appro- 
priate for an Scd galaxy near z=0.2 of 0.45 mag (Frei & 
Gunn 1994), N9-2 has /^b(0) = 22.1 well within the range 
of high surface brightness galaxies. N9-2 is just a normal 
HSB spiral galaxy, seen at moderate redshift, and so we 
remove it from the rest of the analysis. 

The galaxies NlO-4 and C5-5 (Figure 3e) both have 
Mb(0) = 24 and are among the faintest in our sample. The 
S/N of the observations are sufficient to determine red- 
shifts for these two galaxies, but insufficient for a more de- 
tailed analysis. NlO-4 is particularly interesting, however, 
because of its companions. There were two other objects 
which serendiptously fell onto the slit during the observa- 
tions. We refer to these as N10-4a and N10-4b. N10-4a 
(ll''58"52^5, 20°58'23" J2000) lies 4h-i kpc away (12") 
from the center of NlO-4. We detect emission lines of Ha, 
H^, and [O III]AA4959,5007 in the companion, shifted 65 
km s~^ blueward (cz/je;=7382 km s~^). We do not detect 
anything in the broadband setup image at this position, 
making it difficult to determine whether this is a sepa- 
rate system, or a knot of emission in the outer disk of 
NlO-4. Additionally, there is another companion 16". 5 in 
the opposite direction along the slit (N10-4b, ll'*58™50^8, 
20°58'33".5 J2000). This object is also not seen in the 
setup image, or in images from the digitized sky survey. 
We do not detect any continuum from this object and only 
one emission line is seen. This line is detected in all four 
observations of NlO-4. The spectrum is presented in Fig- 
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Fig. 3b. — Spectra for LSB galaxies with emission lines, ordered by the strength of [N IIJA6584 relative to Ha. 



ure 4. The observed wavelength of this emission line is 
5738.9 A, which translates to 5600 A in the restframe of 
NlO-4. This wavelength does not match any commonly 
observed galaxy emission lines. Instead, we suggest that 
this object is a background galaxy, and the single line ob- 
served at 5738.9 A corresponds to [O II]A3727 seen at a 
redshift of 0.54. At this redshift, our restframe wavelength 
coverage is 2790 A to 4700 A, where we do not expect to 
see any other strong emission lines. 

The SED of C5-5 (Figure 3e) has a bump in it, centered 
at 5600 A (5670 A in the observed frame) which we do not 
have an explanation for. The bright spiral galaxy UGC 



4416 fills the edge of the slit, and some of its light might 
be scattering in. There is also a spatial offset between the 
emission line peak and the peak of the continuum for this 
galaxy, though the emission is still consistent with coming 
from C5-5. Due to the low S/N of these observations, we 
do not include either NlO-4 or C5-5 in the remainder of 
the analysis. 

4. EMISSION LINE SOURCE: H II REGIONS VS. AGN 

We use the diagnostic diagrams of Veilleux & Oster- 
brock (1987) to determine whether the emission lines 
in the observed LSB galaxies come from active galac- 
tic nuclei (AGN) or H II regions. Figure 6a shows the 
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Fig. 3c. — Spectra for LSB galaxies with emission lines, ordered by the strength of [N II]A6584 relative to Ha. 



[N II]A6584/Ha vs. [O IIl]A5007/H/3 diagnostic. We over- side of the diagram. They also he in the same area of the 
plot the theoretical separation line of Kewley et al. (2001a) diagram as the HSB galaxies from the Nearby Field Galaxy 
on the diagram. The LSB galaxies all clearly lie in the H II Survey (NFGS; Jansen et al. 2000). Figure 6b shows the 
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[S II]A6717+A6731/Hq! vs. [O III]A5007/H/? emission line 
diagnostic. The interpretation here is less clear. The LSB 
galaxies all lie on the separatrix, and there is considerably 
more scatter amongst the NFGS galaxies as well. Kcwlcy 
et al. (2001b) note that the [S II]A6717+A6731/Ha ratio is 
sensitive to the local density, and therefore might be a less 
rehable diagnostic than the [N II]A6584/Hq; ratio. NlO-2 
shifts dramatically between the two diagnostics, relative 
to the rest of the LSB galaxies. This is an instrumental 
effect, as the redshifted [S II] lines fall just at the red end 
of our spectral coverage, where the flux calibration is less 
secure. All the galaxies with emission lines arc consistent 
with star formation as the ionizing source, rather than 
AGN. 

5. MEAN AGE AND METALLICITY OF THE STARS 

The degeneracy between the effects of age and mctallic- 
ity on broadband optical colors and individual spectral fea- 
tures has been well documented (O'Connell 1976; Worthey 
1994; Rose 1994). This degeneracy can be broken by us- 
ing combinations of two or more absorption features, with 
different sensitivity to age and metallicity (Worthey ct al. 
1994; J0rgensen 1997; Vazdekis 1999; Kuntschner et al. 
2001). We use the line index H/3 which is more age sensi- 
tive, in conjunction with the more metal sensitive indices 
Mgb and < Fc >. In Figure 7 we plot grids of age and 
metallicity for single stellar populations derived from the 
V2000 models, on the axes of H/3 vs. Mgb and H/J vs. 
<Fe>. Within this grid, we determine the mean luminos- 



ity weighted age and metallicity for the stellar populations 
of the four absorption line dominated LSB galaxies, Cl- 
2, Cl-4, N3-1, and Ul-8. We also plot the two elliptical 
galaxies which we observed, as well as the locus of points 
covered by the Coma cluster elliptical and SO galaxies from 
J0rgensen (1999). 

The three LSB galaxies Cl-2, Cl-4 and Ul-8 all lie to- 
gether on the plots, in the low metallicity regime. They 
all have mean ages around 5.5 Gyr. Cl-4 has near so- 
lar metallicity, while Cl-2 and Ul-8 have about half solar 
metallicity. 

In contrast, N3-1 has a much higher metallicity, appear- 
ing at 1.5 times solar metallicity in the 11/3 vs. < Fe > 
plot, and even higher in the H/3 vs. Mgb plot. This is 
the regime shared by luminous elliptical and SO galax- 
ies. Furthermore, note how N3-1 as well as the ellipti- 
cal and SO galaxies shift with respect to the model grid 
between the two plots. This is the signature of super- 
solar [Mg/Fe] abundance ratios. This enhancement effect 
is well known for luminous galaxies, both spirals and ellip- 
ticals, and the strength of the enhancement is correlated 
with galaxy mass and velocity dispersion (Worthey 1998; 
J0rgensen 1999; Kuntschner et al. 2001). The other three 
LSB galaxies do not move relative to the model grid be- 
tween the two plots; these galaxies all have solar [Mg/Fe] 
abundance ratios. 

Figure 8 more clearly illustrates the [Mg/Fc] issue. This 
plot shows the < Fe > index plotted against the Mgb in- 
dex. Since these two indices have similar metallicity and 
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age sensitivity, the model grid for differing ages and metal- 
licities is not well separated in this plane. Galaxies with 
solar [Mg/Fe] should fall onto the narrow region covered by 
the model grid, as indeed most of the LSB galaxies do. On 
this Figure we also show the LSB galaxies with emission 
lines which have reliable Mgb and < Fe > measurements. 
The bulk of the elliptical and SO galaxies as well as N3-1 
fall to the right of the model grid, a failure of the models 
to account for the [Mg/Fe] appropriate to these galaxies. 

There is another caveat with these models, especially 
important when dealing with low metallicity stellar pop- 
ulations. These models do not include blue horizontal 
branch or post-AGB stars, which will have an effect on 
measured Balmer line strengths (Maraston & Thomas 
2000). However, at the lowest metallicities found here (a 
tenth solar to a fifth solar), galactic globular clusters do 
not show a blue horizontal branch, and thus we do not ex- 
pect contamination by these blue stars to be a significant 
effect. Furthermore, the effect of these stars on integrated 
spectra is to increase the strength of the H/? absorption, 
leading to younger apparent ages. The main conclusion 
of this section is that some LSB galaxies do in fact have 
old mean ages; a correction for these blue stars would only 
strengthen our conclusion. 

6. GAS PHASE METALLICITY AND AGE INDICATORS 

Much as we used the stellar absorption indices to study 
the older stellar populations, we would like to use emission 



line strengths and ratios to characterize the metallicity of 
the gas and the properties of the young stellar popula- 
tions. The equivalent width of Ha, EW(Hq:), compares 
the flux from hot young stars, capable of ionizing nebu- 
lae, to the red flux from the stellar population, which is 
dominated by long-lived lower mass stars for a standard 
(e.g. Salpeter 1955) initial mass function (IMF). Conse- 
quently, the EW(Ha) can be used as a mean age indicator 
for the integrated stellar population. The conversion be- 
tween EW(Ha) and the Scale (1986) birthrate parameter, 
defined as 6 = <sFR>past ' '® sensitive to both the form of 
the star formation history (SFH) and the IMF, and hence 
is quite model dependant. A single burst of star forma- 
tion will only produce Ha emission for ~ 30 Myr after 
which the massive stars (M^ > 15 M©) which produce 
the ionizing flux have all died off. Models with continuous 
star formation and multiple burst models that produce an 
underlying old population plus a new burst will both pro- 
duce Ha emission at later times. Kennicutt, Tamblyn, & 
Congdon (1994) use stellar population and nebular mod- 
els to derive a relation between EW(Ha) and h appropri- 
ate for galaxies with a constant or exponentially decreas- 
ing SFR. The relation they derive, based on a Salpeter 
(1955) IMF with upper and lower mass limits of 100 and 
0.1 Mq, a 10 Gyr age and varying exponential decay times, 
shows a monotonic increase of h with EW(Ha), and has 
6=1 at EW(Ha)=64 A. van den Hoek et al. (2000) found 
that chemo-evolutionary models with an exponentially de- 
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creasing SFR were sufficient to describe the LSB galaxies 
in their sample with (B— V)>0.4 but that additional star 
formation bursts were needed to explain bluer galaxies. 
The measured EW(Hq;) for an old-plus-burst stellar pop- 
ulation will of course depend on both the age of the burst, 
and the stellar mass fraction that it embodies. In Figures 
10 and 12 we show the b calibration of Kennicut et al. 
(1994) on the right vertical axis, but give the caveat that 
this calibration is inappropriate for an old-plus-burst SFH. 

Traditional methods of studying gas phase chemical 
abundances in H II regions require measurements of sev- 
eral H or He recombination lines, along with collision- 
ally excited lines from multiple ionization states of a 
heavier metallic species. Ostcrbrock (1989) details the 
method for studying individual H II regions. Kobul- 
nicky et al. (1999) extend this work and look at the ef- 
fect on the derived abundances of smearing together the 
light from multiple H II regions in a single galaxy into 
a single spectroscopic aperture, using abundance analy- 
sis methods such as the semi-empirical strong-line R23^ 
method (Pagel et al. 1979; Osterbrock 1989; McGaugh 
1991, 1994). The wavelength range covered by our spec- 
tra rarely includes the [O II]A3727 A fine used in the R23 
method. Instead, we choose to use the metallicity indi- 
cator [N II]A6584/Ha (hereafter [N II] /Ha) proposed by 
Storchi-Bergmann, Calzetti, & Kinney (1994). This mea- 
surement has the advantage of being monotonic (as op- 

4 log(R23)= log('°"l^^^^^+t°'"l^^^^^^-^°°^) 



posed to the double valued nature of the R23 method) and 
is essentially reddening independant by virtue of the small 
wavelength split between the two lines. Empirical calibra- 
tions between [N II] /Ha and log(0/H) (determined us- 
ing the R23 method) have been derived from observations 
of individual H II regions (e.g. van Zee et al. 1998) and 
for galaxies (Denicolo, Terlevich, & Terlevich 2002). The 
drawback to using [N II] /Ha is that this ratio is sensitive 
to the ionization state, which will vary between H II re- 
gions and the diffuse ISM, and is difficult to determine for 
globally averaged spectra. Kewley & Dopita (2002) pro- 
vide a theoretical calibration of [N II] /Ha vs. metallicity 
for starburst galaxies, including the effects of the ioniza- 
tion parameter. McGaugh (1994) find for H II regions 
in LSB galaxies that the ionization parameter log<U> 
varies between —2 and —3.5, while Denicolo et al. find 
that log<U>=— 2.0 is most representative for their sam- 
ple of galaxies. However, a change in log<U> of -1-0.5 dex 
or —1.0 dex leads to an uncertainty in the metallicity zero- 
point of ±0.3 dex. There was no correlation seen between 
metallicity and log<U> in the McGaugh dataset. In the 
plots that follow, we label one axis with the Denicolo et al. 
empirical calibration of [N II] /Ha, but we give the caveat 
that the zeropoint of this calibration may be high by ~0.3 
dex and that any systematic changes of the ionization pa- 
rameter with other parameters of these galaxies such as 
age may lead to a misinterpretation of the [N II] /Ha rank- 
ing as a strict metallicity scale. 
We compare the gas phase metallicity indicator [N II] /Ha 



21 



-20 20 




Distance From Galaxy Center (arcsec) 



Fig. 5. — The galaxy Ul-4: (a) image of Ul-4 taJjen with the 
HET during set-up for the spectroscopic observation. The object 
at the lower right is a bright star. The three arrows correspond to 
the regions of the slit where emission peaks arc seen, (b) spatial cut 
across the 2D spectrum. Solid line - flux of the Ha line (average over 
6569. oA - 6565.9A in the Ul-4 restframe); Dashed line - continuum 
flux blueward of Ha (average over 6495.3A - 6540.8A in the Ul-4 
restframe). The Ha emission is more spatially extended than the 
stellar continuum light, and three peaks arc scon in the emission 
(spatial) profile, at -12", -5" and +5" from the galaxy center. We 
don't see the emission regions in (a), which is a broadband image. 
The spatial profile of the star in (b) shows the image PSF for this 
observation (the wavelength regions in the restframe of the star are 
6633.7A - 6640.7A and 6569.3A - 6615A). 



with the stellar absorption indices Mgb and < Fe > in Fig- 
ure 9. There are clear correlations with each index. A 
Kendall's r test gives a probability of 1.6% that [N II]/Ha 
and Mgb are not correlated, and a probability of 0.08% 
that [N II] /Ha and <Fe> are not correlated. The best fit 
relations are: 

[NII]/Hq! = -0.11(±0.13) + 0.26(±0.08) • <Fe> (1) 
[NII]/Ha = -0.01(±0.09) + 0.16(±0.05) ■ Mgb (2) 

It is perhaps surprising that the correlations are so strong, 
given the age-metallicity degeneracy of the stellar indices 
and the variation of the [N II] /Ha calibration with ion- 
ization parameter. Panels c and d in Figure 9 show 
the age indicator EW(Ha) plotted against the absorption 
line indices. At low Mgb strength, we sec a wide range 
of EW(Ha) strengths, which diminishes at higher Mgb 
strengths. In models with continuous, exponentially de- 
creasing SFR, the metallicity increases as the gas fraction 
decreases. Simultaneously, the SFR drops and the frac- 
tional galaxy mass made up by low mass stars increases. 
This will result in a decrease of EW(Hq;) and simultaneous 
increase of Mgb over time, leading to an upper envelope 
in the EW(Ha) vs Mgb plane, much as we see here. For 
"open box" chemical evolution models (models in which 



enriched gas is allowed to escape due to stellar winds and 
supernovae), the low EW(Ha) and low Mgb region will 
get filled in by galaxies with low total mass, where the 
star formation rate drops once the gas is expelled, and be- 
fore metallicity builds up. In any case, once star formation 
stops completely, the EW(Ha) will drop to zero, and the 
Mgb index will slowly strengthen as the stellar population 
ages (at constant metallicity). This picture holds equally 
well for the < Fe > strength as for Mgb. A quantitative 
analysis of where models with continuous star formation 
would lie on the EW(Ha) vs. metal absorption index plots 
requires models that contain both nebular emission spec- 
tra, and stellar absorption line strengths. Unfortunately, 
there are not yet any such models in the literature. Within 
the picture of continuous star formation, we can also give 
a qualitative explanation for the strong correlation seen 
between [N II] /Ha and Mgb or <Fe>. [N II] /Ha reflects 
the current metallicity of the gas. The ionization param- 
eter is a small source of scatter, but is not (apparently) 
significant. The Mgb and <Fe> strengths arc responsive 
to both the metallicity and age of the stellar population, 
but since these are linked by the gas consumption and en- 
richment history, a correlation exists between the current 
metallicity and the age and enrichment history of the stel- 
lar population. 

An alternative explanation can be given by a model with 
an old stellar population underlying a young burst of star 
formation. In this case, the mass fraction of the burst 
and metallicity of the old population will determine the 
evolutionary track of a galaxy through the EW(Ha) vs. 
Mgb (or <Fc>) plane. The mass fraction determines how 
much blue continuum light from the young stellar popu- 
lation there will be to veil the strength of the underlying 
metal line absorption. The Ha flux is driven by very short- 
lived stars, while the continuum veiling effect evolves over 
longer timescalcs (~ 0.5 Gyr vs. ~ 30 Myr for the Ha). 
Consequently, the tracks will move vertically down the di- 
agram, and then to the right at late times. The track will 
end at the absorption line strength of the underlying pop- 
ulation for small bursts, and the offset between where the 
track descends vertically and where it ends depends on the 
mass fraction of the burst. To illustrate this, we use the 
Starburst99 models (Leitherer et al. 1999) for the evolu- 
tion of an instantaneous burst population, along with the 
V2000 stellar population models for the underlying pop- 
ulation. In Figure 9(c-d) we plot the data for the LSB 
galaxies together with model tracks for galaxies with a 5 
Gyr old, solar metallicity underlying population, and a 
10%, 2%, 1%, and 0.1% (by mass) instantaneous burst. 
It is clear that by varying the age and metallicity of the 
underlying population, or by varying the burst mass frac- 
tion, it is possible to flll in every part of the plot. While 
we do not find LSB galaxies in all parts of the plot, we 
have a small sample and do not strongly constrain against 
the possibility of galaxies in the high EW(Ha) and high 
Mgb (or <Fe>) region of the diagrams. 

It is possible to fully explain the EW(Ha), Mgb, and 
< Fe > data with an old (> 5 Gyr) solar metallicity pop- 
ulation and varying burst mass fractions. However, with 
this model the [N II] /Ha vs. Mgb correlation becomes dif- 
ficult to explain. The [N II] /Ha value should reflect the 
metallicity of the burst, which is presumably greater than 
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Fig. 6.- Ionization diagnostic diagrams: large boxes — LSB galaxies; small boxes — galaxies in the NFGS sample (Jansen et al. 2000); lines 
— theoretical classification line and ±0.1 dcx error range separating AGN from starburst (H II) galaxies, from Kewlcy ct al. (2001a). The 
LSB galaxies are all well within the H II region for the [N II] /Ha vs [O IIIJ/H/3 diagnostic. There is a shift in the data with respect to the 
theoretical separatrix in the [S II] /Ha vs [O IIIJ/H/3 diagnostic. However, the LSB galaxies arc all still consistent with starburst like spectra. 
The shift of NlO-2 between (a) and (b) is a result of poor flux calibration at the extreme red end of the spectrum, which affects the [S II] 
measurements. The arrow shows the uncertainty in [S II]/Ha for NlO-2 due to this systematic effect. 



or equal to tlic metallicity of tlie underlying stellar popula- 
tion. If the underlying stellar population is always metal 
rich and the burst mass fractions were random, pushing 
the absorption line indices to lower levels, there would not 
be a correlation between [N II] /Ha and the absorption in- 
dices. Consequently, it appears that old-plus-burst mod- 
els are only viable if the old populations have a range of 
metallicities, and if any one burst constitutes only a small 
fraction of the stellar mass (less than ~ 0.3%) such that 
the absorption indices are not perturbed from that of the 
underlying population. 

Our understanding of these correlations is hampered 
by the small size of our sample, and discontinuity be- 
tween models which predict nebular emission strengths 
from young stellar populations, and models which predict 
the stellar metal absorption line indices for older, cooler 
populations. Ideally, we could gleen a better understand- 
ing of these correlations from a larger sample of LSB or 
HSB galaxies. However, we have been unable to find pub- 
lished studies of HSB galaxies which include data on both 
the stellar absorption indices and emission line strengths. 
The NFGS and Sloan Digital Sky Survey (though the lat- 
ter is prone to aperture effects) represent two datascts 
upon which it would be enlightening to make this study. 

We cannot measure gas phase abundances for the four 
absorption line dominated galaxies, but we use the correla- 



tions in Equations (1) and (2) to assign values of [N II]/Ha 
to these galaxies. We do this primarily so that we can 
include these four galaxies on the remainder of the diag- 
nostic diagrams, where we use [N II] /Ha as a metallicity 
indicator. We interpolate the relations get values for CI- 
2 and Ul-8, and extrapolate to get values for Cl-4 and 
N3-1. The values determined from each of the two cor- 
relations are averaged, to reduce the uncertainties. The 
assigned [N II] /Ha values are as follows: 0.47±0.15 for 
Cl-2; 0.47±0.14 for Ul-8; 0.55±0.16 for CI-4; 0.71±0.19 
for N3-1. These four galaxies are shown as open circles on 
Figure 9. The absorption line dominated galaxies seem to 
lie at one end of a continuous sequence of spectral proper- 
ties, as illustrated by Figs. 3a - 3d. This suggests that the 
extrapolation of these correlations is reasonable. However, 
we do note that stellar populations at constant metallicity 
will evolve to higher Mgb and < Fe > with age once the 
star formation stops. Though we do not see a large age 
spread among the four absorption dominated systems, age 
will affect the extrapolation of these relations to higher 
index strengths. By including both the absorption and 
emission line galaxies in Figures 10 and 12 we attempt to 
show these galaxies potentially have very similar evolution- 
ary histories (though they are perhaps at different stages 
of that evolution). The emission associated with a small 
burst of star formation is very short-lived, and though the 
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Fig. 7. — H/9 absorption line index vs. the Mgb index (left) and the < Fe > index (right). Boxes: LSB galaxies; Large triangles: our data 
of HSB ellipticals; small triangles: Coma cluster elliptical and SO galaxies from j0rgcnscn (1999). The model grid shows lines of constant 
age (dotted) and constant mctallicity (long dash), from V2000 as labelled. The same age and mctallicity values are used for the grid in both 
panels. The four LSB galaxies plotted are (in order of increasing Mgb, < Fe > strength) Ul-8, Cl-2, Cl-4, and N3-1. Note that the points 
for Cl-2 and Ul-8 lie virtually on top of eachother. The Coma cluster galaxies, N.3-1, and NGC 3872, shift relative to the model grid in the 
two plots. This is caused by the well known super-solar [Mg/Fe] ratio in luminous elliptical and SO galaxies, which is not reproduced by the 
models. 



emission line flux can be significant, the associated effect 
on the chemical evolution and continuum luminosity could 
be slight, and so understanding these two "classes" of LSB 
galaxies, absorption and emission line galaxies, is necessar- 
ily intertwined. 

In Figure 10 we plot [N II] /Ha against the EW(Hq) 
for the LSB galaxies. The mctallicity and birthrate cali- 
brations discussed at the beginning of this section are in- 
cluded on the axes opposing the measured parameters. We 
also plot the HSB galaxies from the NFGS (Jansen et al. 
2000), which satisfy the condition that EW(Ha)>9 A (i.e. 
those which have a small measurement uncertainty for the 
[N II] /Ha line ratio). The LSB galaxies fill the same re- 
gion of the [N II] /Ha vs. EW(Ha) parameter space, sug- 
gesting a similar star formation and chemical enrichment 
history to that of HSB galaxies. At low metallicities there 
is a wide range of mean ages, however the relative current 
SFR decreases as metallicity increases. This makes sense 
for continuous star formation, as the bulk of the gas in 
these galaxies must have been previously processed into 
stars to raise the mean metallicity, and consequently there 
is less gas around to make new stars. 

7. DUST CONTENT 

The amount of internal reddening in LSB galaxies is 
not well known, but of great importance to attempts at 



modelling their formation and evolution. In Figure 11a 
we plot the Balmer ratio, Hq/H/3, versus the metallicity 
indicator [N II] /Ha. The Balmcr ratio can be used to 
derive the internal reddening of the emission line regions. 
Using the Cardelli extinction law (Cardelli et al. 1989), 
with Rv=3.1, we derive the following relation: 



E(B-Ie) = 2.375 ■E(B-V) = -2.65 (3) 

This relation is used to place an E(B— Ic) scale on the right 
axis of Fig.ll. We use the NFGS (Jansen et al. 2000) as 
a comparison sample of HSB emission line galaxies. The 
NFGS galaxies were not observed in the I-band. We con- 
vert their (B— Rc) colors to (B— Ic) colors using the re- 
lation (B-Ic) = 1.34 • (B-Rc) + 0.18, derived from the 
mean galaxy colors of Fukugita, Shimasaku, & Ichikawa 
(1995). Central surface brightnesses were not derived for 
the NFGS galaxies. In the LSB data, we have not cor- 
rcctcxl the measured Ha fluxes and equivalent widths for 
underlying stellar absorption. This correction will move 
points upwards and to the left in the Ha/H/3 vs. [N II] /Ha 
plane. While the uncorrected data show a trend for LSB 
galaxies to have lower extinction than HSB galaxies at the 
same metallicity, an Ha equivalent width correction of 1.5 
A would move most of the LSB galaxies onto the locus of 
HSB galaxy points. 
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Fig. 8. — The < Fe > index vs. the Mgb index. Filled boxes: 
LSB galaxies with no H/3 emission; Open boxes: LSB galaxies with 
measured H/3 emission; Large triangles: our data of HSB ellipti- 
cals; small triangles: Coma cluster elliptical and SO galaxies from 
j0rgensen (1999). The same models used for the grid in Figure 7 
are used here. The Mgb and < Fe > indices have similar sensitiv- 
ity to age and metallicity, so the model grid is very narrow. N3-1, 
like many of elliptical galaxies, lies off the model grid, suggesting 
super-solar [Mg/Fe] abundance. See the text for details. 



McGaugh (1994), in his study of LSB galaxy H II re- 
gions, found evidence for dust in some galaxies, with in- 
ternal extinctions as high as E(B— V) ~ 0.5 (corresponding 
to E(B— Ic) ^ 1.3). The galaxies in that study were sim- 
ilar to the lower metaUicity galaxies ([N II] /Ha < 0.2) in 
Figure 11a for which we find a similar range of redden- 
ings. His study did not include any of the red or giant 
LSB galaxies for which we find generally larger reddening 
values. 

The dust distribution in galaxies tends to be clumped, 
and these reddening values apply only to the regions where 
the emission lines are produced. The emission line regions 
and light from the stellar disk are not spatially resolved 
in our observations and both contribute to the extracted 
spectra. Consequently, we cannot constrain the redden- 
ing aff'ecting the general stellar disk, which dominates the 
continuum luminosity of the galaxy, and thus we do not 
make a correction to the global photometry for internal 
extinction. 

Figure 11a shows a trend for the Balmer ratio (and pre- 
sumed reddening) to increase with increasing gas-phase 
metallicity. In Figure 1 1 (b-d) we plot the Balmer ratio ver- 
sus photometric parameters dominated by the stellar disk. 
There is no correlation at all between Ha/Hfi and central 
surface brightness over the range 22< ^B(0)cor < 25, nor 
is there a correlation between Ha/H/3 and galaxy incli- 



nation. The correlation between (B— lc)tot color and the 
Balmer ratio is shown for the NFGS galaxies in Figure 
lid. The LSB galaxies fall onto the same locus of points, 
especially after applying a correction for underlying Ha 
absorption of 1.5 A. The LSB galaxies appear to lie along 
the reddening vector in Figure lid. This might imply 
that the dust is not local to the emission region, but is 
instead mixed throughout the galaxy. Alternatively, it is 
possible that the colors are not caused by dust redden- 
ing of the stellar light, and the relation between Ha/H/3 
and (B— Ic)tot is not the fundamental one. Rather, the 
amount of dust increases with increasing metallicity (Fig. 
11a), and the stellar populations also become redder with 
increasing metallicity (cf. the tight correlation between 
[N II]/Ha and (B— lc)tot shown on Fig. 12c, and the corre- 
sponding discussion in Section 8), which in turn results in 
the correlation seen here between Ha/H/3 and (B— lc)tot- 
If the colors are driven by dust reddening over the whole 
galaxy, then we would expect to also see a correlation be- 
tween galaxy inclination and reddening, with higher red- 
dening measurements for galaxies closer to edge-on. We do 
not detect any such correlation, and thus conclude that the 
correlation between Ha/H/3 and (B— lc)tot is secondary, 
and both parameters are driven primarily by the metallic- 
ity. 

Regardless of the interpretation of the galaxy colors, it is 
clear that dust is present in abundance in some LSB galax- 
ies, particularly the red and high metallicity ones, but the 
absence of a correlation between dust and /UB(0)cor sug- 
gests that dust is not the cause of the low surface bright- 
ness for LSB galaxies of any color. 

The relative emission line strengths presented in Table 
7 have not been corrected for internal extinction within 
these galaxies. The line ratios we use in our analysis (e.g. 
[N II] /Ha) all have narrow wavelength spreads, and are 
not significantly affected by dust in the amount we de- 
tect. The uncertainty of a correction for intrinsic redden- 
ing, however, would be significantly larger than the cor- 
rection itself (which is negligible), and so we do not apply 
any correction. 

8. EMISSION LINE PROPERTIES VERSUS GLOBAL 

PHOTOMETRIC PROPERTIES 

In Figure 12, we show [N II]/Ha and EW(Ha) ver- 
sus central surface brightness, absolute magnitude, and 
(B— lc)tot- We plot the absorption line galaxies in this 
figure as well, using the [N II] /Ha values derived and dis- 
cussed in Section 6. However, these four galaxies are never 
included when we measure correlations between the pa- 
rameters plotted in Figure 12. 

In Figures 12c and 12f, the outher P9-4 at (B-lc)=2.33 
(OBSC197) is purportedly the reddest object in the sam- 
ple, despite its high EW(Ha) and metal poor nature. The 
continuum slope in our HET spectrum of this object is 
very flat. In comparison, the galaxies at (B— lc)=1.3 have 
bluer continua, but those at (B— lc)=1.8 have redder con- 
tinua. The continuum slope of NlO-2 (B— lc)=1.48 is also 
very flat. We suggest that there are errors in the pub- 
lished photometry for P9-4. We label P9-4 in all panels of 
Figure 12 and we do not include the data for P9-4 in the 
determination of any of the correlations discussed below. 

Figures 12a and 12d show that there is a broad range of 
both [N II] /Ha and EW(Ha) values for the whole range 
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Fig. 9. — Comparison of gas phase indicators and stellar absorption line indiees. The gas phase [N II]/Ha ratio plotted against (a) the 
stellar Mgb line index and (b) the <Fc> line index. The best fitting relations are shown as solid lines, (c) The EW(Htt) vs. the Mgb index, 
(d) EW(Ha) vs. the < Fe > index. Boxes: LSB galaxies with measured emission lines; circles: LSB galaxies with no measured emission, 
see Section 6 for details. The lines in panels (c) and (d) are evolutionary tracks for an old-plus-burst stellar population, with a 5 Gyr old 
population, and bursts which comprise (from left to right) of fO% (dotted), 2% (short-dash), 1% (long-dash) and 0.1% (dot-dash) of the total 
stellar mass. The kink in these evolutionary tracks occurs at ~ 7.5 Myr after the burst, and the EW(Hq) drops below 1 A within 20 Myr 
of the burst for all the tracks. The tracks all end at Mgb = 3.6 and < Fe > = 2.7, EW(Ha) ~ 0, which are the values of the underlying 
population. 



of surface brightnesses studied. A Kendall's r test gives 
a probability of 54% that there is no correlation between 
EW(Ha) and )UB(0)cor- There is a probability of 31% that 
there is no correlation between [N II] /Hq and /LiB(0)cor- In 
contrast, Bell & de Jong (2000) found correlations both 
between mean age and K-band central surface brightness, 
and between mean metallicity and K-band central surface 
brightness. Their modelling involved an asstimption of 
exponentially decreasing star formation rates all begin- 
ning at the same time in the past (12 Gyr), with differ- 
ent exponential timescales, r, which thus produced stel- 



lar populations with different luminosity-weighted mean 
ages at the current epoch. They varied r and metallic- 
ity in their models to produce a grid in the (Rc— K) vs. 
(B— Ic) plane, which was then used to interpret their data. 
These models have a very similar SFH to the models of 
Kcnnicut ct al. (1994), which wc use for the right side 
scale on Figure 12(d-f) and which predict a monotonic 
relationship between EW(Ha) and t, and consequently 
mean age. However, as mentioned above, we do not see a 
correlation between EW(Hq;) and B-band central surface 
brightness. Bell & de Jong make an unweighted least- 
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Fig. 10. - The metallicity indicator [N II]/Ha plotted against 
EW(Hq), a measure of the birthrate parameter. Large boxes: LSB 
galaxies with measured emission Unes; circles: LSB galaxies with 
no measurable emission lines whose positions are based on the ex- 
trapolation described in Section 6; small boxes: HSB galaxies from 
the NFGS (Jansen et al. 2000). The right axis shows the [N II] /Ha 
metallicity calibration of Denicolo et al. (2002), appropriate for H II 
regions. The top axis shows the EW(Ha) calibration of the Scalo 
birthrate parameter for one star formation model detailed in Ken- 
nicut et al. (1994). See text for details. 



squares fit to their data which produces a change in mean 
age of 0.84 ± 0.08 Gyr for a change of two magnitudes 
in central surface brightness, equivalent to the range cov- 
ered by our sample, with lower luminosity-weighted mean 
ages occuring at lower central surface brightness. Their 
fit to the metallicity-surface brightness relation produced 
a change in the mean metallicity of 0.26 dex for the same 
two magnitude change in central surface brightness, with 
higher metallicity at higher surface brightness. Our sam- 
ple size, though similar to the LSB galaxy study of Bell et 
al. (2000), is a small fraction of the combined HSB-plus- 
LSB galaxy sample of Bell & dc .Jong. It is quite possible 
that our sample selection, which deliberately covers a wide 
range of galaxy colors and surface brightnesses, has led to 
more scatter than would be present in a volume-limited 
sample of galaxies, and that therefore we might not be 
sensitive to these predicted modest changes in mean age 
and metallicity. 

The line ratio [N II] /Ha versus the absolute magnitude 
is shown in Figure 12b. This plot makes clear the point 
that LSB galaxies are not synonomous with dwarf galaxies, 
but cover a range of absolute magnitudes comparable to 
that of HSB galaxies. The emission line LSB galaxies show 
a correlation between [N II] /Ha and absolute magnitude. 
A Kendall's r test gives a probability of 5% that no cor- 
relation is present. The relation for LSB galaxies matches 
that of the NFGS galaxies. The LSB galaxies without 



emission also follow the correlation, but form a ridge on 
the faint, high metallicity edge of the galaxy distribution. 
This may be reflecting a higher B-band mass-to-light ra- 
tio in these older objects, if the [N II]/Ha-MBtot relation 
is a result of a more fundamental metallicity-mass rela- 
tion. Alternatively, it could reflect a systematic error in 
the interpolation/extrapolation of the [N Il]/Ha-Mgb and 
[N II]/Hq:-< Fe > relations used for the non-star-forming 
galaxies. 

The distribution of the LSB galaxies is a good match to 
that of the NFGS galaxies in Figure 12e. We do not detect 
any correlation between EW(Ha) and Metot • A Kendall's 
T test gives a probability of 63% that there is no correlation 
for the emission line LSB galaxies. A rough scaling of 
the absolute SFR can be determined by multiplying the 
EW(Ha) by the total R-band luminosity, which includes 
the continuum near Ha. Within our sample, the galaxies 
with the largest (EW(Ha) ■ Lb) arc the ones at the higher 
surface brightness end of our range, with /XB(0)cor ^ 23.0 
mag/n". Galaxies with lower total SFR will be found in 
the lower right of the EW(Hq) vs. Metot diagram. Our 
sample selection is neither complete nor representative, so 
we cannot comment on how the absolute star-formation 
rates of LSB and HSB galaxies compare, van den Hoek 
et al. (2000) found that the mean present-day SFR for a 
sample of late-type LSB galaxies was about one-tenth the 
SFR of similar HSB galaxies. 

The relations between [N II] /Ha and EW(Ha) with 
(B— Ic) color, shown in Figs. 12c and 12f, are compara- 
ble in strength to those between [N II] /Ha and the stellar 
absorption indices Mgb and < Fe >. A Kendall's t test 
applied to the emission line LSB galaxies gives a probabil- 
ity of 2% that there is no correlation between [N II] /Ha 
and (B— Ic)tot- Similarly, there is a probability of only 
0.7% that EW(Ha) and (B-Ic)tot are not correlated for 
the emission line LSB galaxies. (B— Ic) color is sensitive 
to changes in age, metallicity, and dust reddening, and 
all three may play some role in explaining the correla- 
tions. If we use the Denicolo et al. (2002) calibration for 
[N II] /Ha, the resulting metallicity range spanned by the 
data is less than one dex. At any constant age less than 5 
Gyr, the V2000 stellar population models predict a max- 
imum (B— Ic) color range of only 0.6 mag for a range in 
[Fe/H] from -1.0 dex to -1-0.2 dex. The (B - Ic)tot range 
spanned by our data is more than one magnitude, so vary- 
ing only the metallicity of the stellar population is not a 
sufficient explanation for the colors. Stellar populations 
with mean ages older than 5 Gyr will show a wider color 
range over that range of metallicity, but are not supported 
by the ages determined in Section 5 for the four absorption 
line dominated LSB galaxies. 

Including age in the interpretation of the colors requires 
us to make assumptions about the SFH in these galaxies. 
In Section 6, we discussed three plausible star formation 
scenarios: (1) continuous, but decreasing star formation 
rate with a build-up of metallicity, (2) old-plus-burst mod- 
els with an old stellar population underlying a significant 
(> 1% by mass) burst of star formation, and (3) old-plus- 
burst SFH models where the burst was small (< 0.3% 
by mass) and produced Ha emission, but had little effect 
on the global photometric properties. Scenario (2) was re- 
jected because it would not maintain the strong correlation 
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Fig. 11. — The Balmcr ratio, Ha/H/3, plotted against (a) the metallicity indicator [N II] A6584/Htt, (b) the central surface brightness 
/iB(0)cor, (c) galaxy inclination (edge-on = 90°), and (d) The total (B— Ic) color. Large boxes - LSB galaxies; small points - galaxies from the 
NFGS (Jansen et al. 2000). The right side y-axis lists the E(B— 1^) values derived from equation 3. The vertical arrows show where the LSB 
galaxies would move if we made a correction of l.-S A to the Ha equivalent width to account for underlying stellar absorption. The dashed 
arrow in panel (d) shows the locus of reddening vectors for the LSB galaxies. The reddening vector for the NFGS galaxies (whose colors have 
been transformed from observed (B— Rc) colors) would be ~10% steeper. The uncertainties shown for the Balmcr ratios do not include those 
caused by uncertainties in the absorption correction to H/3. The NFGS data has not been corrected for underlying stellar absorption at Ha 
or H/3. Applying an arbitrary correction of 1.5 A to the equivalent width of both Ha and H/3 does not significantly shift the location of the 
locus of HSB galaxy points. 



between [N II] /Ha and Mgb or <Fe>. Scenario (3) was 
acceptable simply because the blue flux from the young 
population did not significantly veil the strength of the 
underlying metal absorption strength, and the [N II] /Ha - 
metal index correlations were unperturbed. However, this 
assumption also means that the global colors will not be 
significantly changed by the young burst. Consequently, 
for this SFH we must invoke a screen of dust as the expla- 
nation for the large range of (B— Ic) colors. Scenario (1) 
predicts a correlation between age and metallicity, which 
then produces a tight correlation between metallicity and 
the Mgb, < Fe > strengths for the stellar population. By 
invoking a correlated evolution of age and metallicity, we 
avoid adding scatter to parameters affected by both, in- 
cluding color. Furthermore, the range of (B— Ic) colors 
in this picture is at least as great as the range seen in 
our data. The V2000 models for a 1 Gyr population 
with one-fifth solar metallicity predict (B— Ic)=1.57, while 
a 4 Gyr solar metallicity stellar population would have 
(B— Ic)=2.55. Mean ages less than 1 Gyr yield colors bluer 
than (B— Ic)=1.5, and while the zeropoint of the model col- 
ors is poorly calibrated, the range is more accurate. Dust, 



while certainly present in the emission line regions, is not 
necessary to explain the global colors in this scenario. 

9. DISCUSSION 

The previous sections show that LSB galaxies cover a 
wide range of metallicities, and as a class they may contain 
both old and young stellar populations. Dust is present 
in some LSB galaxies, and may affect their global col- 
ors. Star formation occurs continuously in these galaxies, 
but with global rates which have varied exponential decay 
timescales. The metallicities of the galaxies are coupled to 
their ages, such that galaxies with older mean ages have 
reached a higher metallicity. Bursts that occur after long 
periods without star formation must make up a very small 
fraction of the total galaxy mass. 

The detailed study of LSB galaxies is still a young field, 
partly because the observations are very difficult to make. 
The study of LSB galaxies originated with the realization 
that selection effects in imaging surveys had an impact on 
our understanding of the galaxy population, and galaxy 
formation and evolution in general. While we are now 
more conscious of some of the biases, samples chosen for 
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follow-up study are still subject to selection effects. It is 
interesting to look at the studies in the literature in the 
context of our results, and determine what effects selec- 
tion or observational biases have had on both our and the 
published conclusions. 

Previous spectroscopic studies of LSB galaxy chemical 
abundances (McGaugh 1994; Ronnback & Bergvall 1995; 
de Blok & van der Hulst 1998) were done on galaxies se- 
lected to have at least one bright H 11 region. These stud- 
ies concluded that LSB galaxies were all metal-poor, with 
gas phase metal abundance less than one-tenth to one-fifth 
solar. These galaxies were selected from photographic sur- 
veys which were much more sensitive to blue galaxies than 
red galaxies at the surface brightness limit. Furthermore, 
the H II regions selected were the brightest H II regions 
in the galaxy, and lower metallicity H II regions will have 
higher ionization parameters and be brighter than high 
metallicity H II regions with similar numbers and ages of 
stars. Compared to our sample, and based on the implicit 
color criterion imposed by the photographic surveys, these 
studies were limited to having galaxies as blue or bluer 
than the five bluest galaxies in our survey. The five bluest 
galaxies in our study (C5-3, N9-1, Ul-4, C4-1, and Pl-3) 
all have active star formation, and define the low metallic- 
ity end of the sample. Though they do overlap the range 
of surface brightnesses and total absolute magnitudes cov- 
ered by the rest of our sample, they are certainly not rep- 
resentative of the full continuum of LSB galaxy properties 
found for the sample. 

Gerritsen & de Blok (1999) made an N-body simulation 
of the SFH for one of the blue late-type LSB galaxies from 
de Blok & van der Hulst (1998), and found that they re- 
quired both low mass surface density and low metallicity 
to explain the properties of the galaxy. The SFH they 
derived has a slowly decreasing SFR when averaged over 
timescales of order 100 Myr, but has significant (factors 
of 2) deviations from the mean SFR on timescales of or- 
der 20 Myr. It is these deviations from the mean, they 
say, which explains the color range seen for LSB galaxies 
in general. Their models do not predict absorption line 
strengths, but if the colors arc significantly affected by the 
SFR deviations, then the Mgb and < Fe > absorption line 
indices must also be affected. Their simulation cannot be 
used to explain the full range of LSB galaxy parameters 
that we see. First of all, despite some uncertainty in the 
metallicity zcropoints of the [N II] /Ha, Mgb and < Fe > 
indicators, we find LSB galaxies with metallicities near so- 
lar, and certainly at least half-solar. Their models require 
metallicities less then about one-fifth solar to maintain a 
low SFR. Furthermore, the range of color (and we infer 
also Mgb, <Fe>) is decoupled from increases in metallic- 
ity, in conflict with the [N II]/Ha-Mgb, [N II]/HQ-<Fe>, 
and [N Il]/Hfi-(B— Ic) correlations. Determining whether 
the range of colors (and thus Mgb and < Fc>) they predict 
at any given metallicity could be hidden within the scatter 
of these correlations requires more quantitative modelling 
than we have done. The Gerritsen & de Blok models pre- 
dict that ~ 20% of the late-type LSB galaxy population 
will have red colors, (B-V) ^ 1. However, those galax- 
ies should be gas rich and metal-poor, and while the gas 
fractions of the red LSB galaxies in our sample remain un- 
constrained, their metallicities are higher than predicted 



by Gerritsen & de Blok. In addition, they predict that 
these red galaxies will have /iB(0) — 24.5 mag/n", which 
is fainter than the galaxies selected for our sample. This 
surface brightness is well below the completeness limit of 
the OBC97 survey, and no other large surveys for red LSB 
galaxies have yet been published. 

van den Hoek et al. (2000) modeled the SFH of the 
galaxies studied by de Blok & van der Hulst (1998) using a 
galactic chemical and photometric evolution model. They 
were able to explain the colors, magnitudes, and gas phase 
metallicities of most of the blue LSB galaxies with a SFH 
that had a continous, exponentially decreasing SFR over 
a period of 14 Gyr. For the bluest galaxies in the sample 
(with (B— V) < 0.45) an additional burst of star formation 
was needed. While our analysis rejects the possibility for 
significant late-time, i.e. high metallicity, bursts, none of 
the galaxies in our sample are that blue. 

Bell ct al. (2000) used near-infrared and optical colors 
to study the stellar populations of blue, red, and giant 
LSB galaxies. They did not require that the galaxies have 
H II regions, nor that they had previous detections in H I. 
However, their sample did have a higher average surface 
brightness than the spectroscopic studies, so that they 
could get sufficient S/N in the K-band imaging. Their 
sample has a similar central surface brightness distribu- 
tion to ours. Like us, they found that LSB galaxies cover a 
range of metallicities, with some having metallicities near 
solar. They also found that some red LSB galaxies had 
rather old luminosity- weighted mean ages. They only had 
five red LSB galaxies in their survey and they determined 
that those with old ages and solar metallicities were all 
at the bright end of the LSB surface brightness range, 
A*b(0) ~ 22.5 mag/n". We find old, metal rich stellar pop- 
ulations in some red galaxies with Hb{0) — 24.0 mag/D". 
By combining the Bell et al. LSB sample with a sample of 
HSB disk galaxies. Bell & de Jong (2000) found a clear cor- 
relation between mean age and surface brightness, which 
we do not detect. The discrepancies in both cases can be 
explained by small sample sizes, especially of the reddest 
LSB galaxies in the sample. 

Impcy ct al. (2001) obtained low resolution (= 20A), low 
S/N optical spectra for 250 galaxies selected from the APM 
survey of Impey et al. (1996), 93 of which have central 
surface brightnesses fainter than /xb = 22mag/n"(their 
LSB cutoff). They also obtained HI spectra for 145 LSB 
galaxies. Burkholder et al. (2001) analyze these data, and 
measure Balmer decrements (Ha/H/?) and [0/H] (using 
the R23 method) when the optical spectra are of suffi- 
cient quality (43 HSB and 17 LSB galaxies have measured 
[0/H]). They also measure the best-fit TuUy-Fisher rela- 
tion for the HI sample. Their sample selection criteria were 
similar to those of Bell et al. (2000) and ourselves: no re- 
quirements for prior HI detections, or the presence of HII 
regions. Though the quality of their data for individual 
galaxies may not be as good, they have a larger datasct to 
draw conclusions from, and their conclusions match ours: 
there is evidence for old stellar populations in some LSB 
galaxies, and the range of metallicities for LSB galaxies 
shows a large overlap with the range of metallicities seen 
in HSB galaxies. 

A complete picture of the formation and evolution of 
LSB galaxies is still not in hand, as many of the conclu- 
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Fig. 12. — The mctallicity indicator [N II]/Hq and the age indicator EW(Hq) plotted against central surface brightness, /iB(0)cor, absolute 
magnitude, Matot! (B— Ic) color. Large boxes: LSB galaxies with measured emission lines; circles: LSB galaxies with no measured 
emission lines; small pentagons: emission line galaxies in the NFGS (Jansen et al. 2000). On the right side axes are the mctallicity and 
birthrate parameter calibrations described in Section 6. The correlation between mctallicity and absolute magnitude found for HSB emission- 
line galaxies is also seen in LSB galaxies. We do not detect a significant correlation between metallicity and surface brightness in this sample. 
The strongest correlation is between color and [N II]/Ha. See text for comments on the published photometry for P9-4. 



sions of these authors and otirsclves are clearly still subject 
to variance due to small and incomplete samples. Future 
detailed studies should rely on deep, complete, multi-color 
surveys with well determined detection criteria to gain a 
better overall picture of LSB galaxy formation and evolu- 
tion. The constraints on LSB SFH implied by the many 
correlations we detect are only applicable to LSB galaxies 
with 1.0 < (B - Ic) < 2.1 and ;Ub(0) < 24.3mag/n". 

10. SUMMARY & CONCLUSIONS 



We have used the Marcario Low Resolution Spectro- 
graph on the 9.2m Hobby-Eberly Telescope to obtain deep 
integrated spectra of the gaseous and stellar components 
of 19 LSB galaxies, covering a range of colors and surface 
brightnesses. These spectra have sufficient spectral resolu- 
tion and S/N in the continuum to measure both emission 
line ratios and stellar absorption line strengths. 

The spectra of these galaxies qualitatively resemble the 
spectra of HSB galaxies covering the full range of spectra 
seen in galaxies of Hubble types from E to Irr. In most of 
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the parameters we measure, LSB galaxy properties cover 
the same broad range as HSB galaxies. Four galaxies have 
spectra whose features are dominated by old stellar pop- 
ulations, covering a range of abundances from less than 
half solar to twice solar metallicity. None of these four 
have been detected in H I. The remainder of the galaxies 
show evidence for ongoing star formation, in some cases at 
a significant rate compared to their average past SFR. At 
least one massive red LSB galaxy shows evidence for super- 
solar [Mg/Fe] abundance. This is in accord with similar 
findings for HSB galaxies, and suggests that the mecha- 
nism responsible for the super-solar abundance ratios is 
effective in LSB galaxies as well. All of the emission line 
galaxies are consistent with star formation, rather than 
AGN, being the ionizing source. 

There are tight correlations between the stellar absorp- 
tion line indices and the gas phase indicators [N II] /Ha 
and EW(Hq!). These have not been studied before, in LSB 
or HSB galaxies. We use these correlations to argue that 
the SFH for these galaxies must be fairly smooth. Large 
bursts of star formation, or multiple small bursts separated 
by long quiescent periods, do not fit these relations. 

Dust is present in some of these galaxies, with Ay as 
large as 1.8 mags. We are unable to determine, though, 
whether the dust is isolated to the emission line regions, 
or if it is mixed in with the the stellar population. 

The redshift distribution for this galaxy sample is very 
broad, ranging from 3000 to 70000 kni/s. Additionally, we 
find several discrepancies between the published redshifts 
based on Arecibo 305m telescope H I observations and our 
new optical observations. In most cases, the discrepancy 
can be explained by beam confusion affecting the Arecibo 
detection. All of the galaxies fomid by O'Neil ct al. (2000a) 
to deviate significantly from the Tully-Fisher relation, we 
show to suffer from significant redshift errors. 
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